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Abstract 

This  study  attenpts  to  determine  if  a pilot  can  manually  release 
a bomb  while  in  a turn.  The  nonlinear  equations  describing  the 
geometry,  fire  control  law,  aircraft  equations  of  motion,  flight 
controls  and  pilot  model  are  developed.  These  equations  are  linearized 
so  that  a frequency  response  analysis  can  be  conducted  for  perturba- 
tions about  a nominal  trajectory.  The  system  response  is  evaluated 
using  wind  gust  inputs  and  lateral  stick  inputs  and  plotting  the 
resulting  system  perturbations  over  a given  frequency  range.  The 
results  of  this  analysis  indicate  the  pilot  is  capable  of  manually 
releasing  a bomb  while  in  a turn.  His  performance  which  is  deter- 
mined by  the  magnitude  of  the  resulting  perturbations  is  comparable 
to  the  fully  automated  system.  This  roll  task  does  not  create  an 
excessive  workload  for  the  pilot. 


I.  Introduction 


Background 

During  recent  years  there  has  been  a growing  need  for  a change 
in  bomb  delivery  tactics.  This  is  primarily  due  to  enemy  defenses 
being  equipped  with  radar  guided  weapons.  For  example,  the  new 
Russian  anti-aircraft  guns,  which  use  linear  predictors  as  part  of 
their  aiming  systems,  are  believed  to  be  extremely  accurate  ivhen 
used  against  aircraft  in  a wings  level,  low  "g"  maneuver.  Therefore, 
in  a heavily  defended  environment,  it  is  desirable  that  an  aircraft 
be  continually  making  heading  and  altitude  changes  to  complicate  the 
enemies  tracking  problems. 

Current  bombing  tactics  and  equipment  require  an  aircraft  to 
track  the  target.  This  is  a wings  level,  relatively  low  "g"  maneuver 
which  requires  several  seconds  for  an  accurate  release.  This  maneuver 
falls  v/ell  within  the  enemies  capabilities  for  tracking  and  predicting 
the  flight  path  of  the  aircraft  and  greatly  increases  the  chances  of 
the  aircraft  being  destroyed.  This  gives  rise  to  the  need  for  a new 
approach  to  get  the  aircraft  to  a reliable  release  point. 

This  new  approach  is  currently  being  studied  by  the  Air  Force 
through  a contract  with  the  General  Electric  Company.  A fire  control 
system  (FCS)  is  being  developed  by  GE  under  the  project  name 
"Firefly".  This  FCS  has  the  capability  of  tracking  a target, 
computing  a release  point,  and  directing  the  autopilot  to  fly  the 
aircraft  to  this  release  point  (Fig.  1).  The  flight  path  from  the 
point  where  the  system  is  engaged  to  the  bomb  release  point  is  primarily 
determined  by  the  number  of  "g's"  which  the  pilot  wishes  to  use.  When 
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the  FCS  is  engaged,  the  system  maintains  this  same  load  factor  and 
computes  the  required  bank  angle  to  take  the  aircraft  to  the  computed 
release  point.  If  the  aircraft  is  fairly  close  to  the  target  for  a 
given  load  factor,  the  bank  angle  will  be  relatively  steep  to  accommodate 
a diving  turn  and  hence,  a short  release  range.  However,  if  the  distance 
to  the  target  is  much  greater  for  the  same  number  of  "g's",  the  bank 
angle  will  be  shallow  to  allow  a climbing  trajectory  which  results  in  a 
greater  release  range.  The  only  restrictions  on  the  system  are: 

(1)  For  a particular  load  factor,  the  aircraft  must  be  close 
enough  to  the  target,  so  that  the  bomb  can  be  "tossed"  into  the 
target  through  a lofted  trajectory. 

(2)  For  a close  in  position,  the  roll  attitude  is  restricted  so 
as  to  ensure  ground  clearance  (Ref  1:5. 1-5.3). 

Purpose 

The  purpose  of  this  study  is  to  determine  if  a pilot  can  manually 
fly  the  aircraft  through  basically  the  same  maneuver  that  would  be 
accomplished  by  the  FCS  coupled  to  the  flight  control  system.  For  this 
analysis,  the  same  target  information  will  be  delivered  to  the  FCS,  however, 
the  output  from  the  fire  control  will  go  to  a heads-up-display  (HUD) 
instead  of  the  flight  controls.  It  will  then  be  the  pilot's  task  to 
determine  an  acceptable  load  factor  and  to  trim  the  aircraft  accordingly. 

He  will  then  be  required  to  follow  the  roll  commands  as  displayed  on  the 
HUD. 

If  this  approach  proves  feasible,  then  the  flight  control  coupler 
can  be  disconnected  from  the  FCS  so  the  pilot  is  in  full  control  of  the 
aircraft  throughout  this  very  critical  maneuver  when  evasive  action  may 
be  required  at  any  time.  This  will  not  only  result  in  a more  responsive 
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aircraft  for  evasive  maneuvers  but  will  eliminate  the  need  for  the  equip- 
ment which  will  tic  the  PCS  into  the  flight  controls. 

Approach 

This  study  is  based  upon  the  system  model  shown  in  figure  2 where 
the  pilot  is  attempting  to  null  the  error  signal  displayed  on  the  HUD. 

The  error  signal  is  a result  of  the  external  inputs  such  as  wind  gusts 
disturbing  the  aircraft  motion  from  some  equilibrium  or  nominal  trajectory. 
As  a result  of  these  disturbances,  perturbation  signals  are  fed  back  into 
the  flight  controls  by  the  stability  augmentation  and  control  augmentation 
systems.  These  disturbances  in  the  flight  controls  along  with  the  pilot 
lateral  stick  inputs  produce  control  surface  deflections  which  in  turn 
affects  the  aircraft  motion.  The  fire  control  law  then  determines  a roll 
rate  command  based  upon  the  aircraft  motion  and  the  position  and  orientat- 
ion of  the  aircraft  relative  to  the  target.  The  roll  rate  command  is  then 
integrated  to  obtain  a bank  angle  command  which  is  displayed  on  the  HUD. 

In  order  to  analyze  this  system,  it  is  first  necessary  to  obtain  the 
nonlinear  differential  and  algebraic  equations  which  represent  the  following 
subsystems: 

(1)  Aircraft 

(2)  Flight  controls 

(3)  Geometry  (Aircraft  position  relative  to  the  target) 

(4)  Fire  control 

(5)  Turbulence  model 

(6)  Pilot  model 

The  geometry  equations  provide  the  nominal  conditions  required  to  define 
the  trajectory  which  will  take  the  aircraft  to  an  accurate  release  point. 
Since  the  nominal  conditions  are  defined  as  an  equilibrium  condition,  it 
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is  appropriate  to  use  the  small  perturbation  technique  to  study  the  system 
response.  Since  frequency  response  is  useful  in  analyzing  the  system 
response,  the  nonlinear  equations  have  been  linearized  by  using  the 
Taylor's  series  expansion  and  retaining  only  the  linear  terms.  The  com- 
puter program  described  in  Appendix  D solves  these  linear  equations  by 
evaluating  the  wind  gust  inputs  or  lateral  stick  inputs  and  plotting  the 
resulting  perturbations  over  a given  frequency  range. 

Knowing  the  lateral  response  of  the  system  and  the  magnitude  of  the 
perturbations  involved,  it  is  possible  to  compare  how  well  the  manual 
system  responds  to  the  fully  automated  system.  Also,  knowing  the  magnitude 
of  the  perturbations  makes  it  possible  to  determine  how  large  a miss 
distance  one  might  expect  due  to  disturbances  such  as  wind  gusts.  Finally, 
a Cooper-Harper  rating  can  be  determined  which  is  indicative  of  the  pilot's 
workload  For  this  roll  task. 

I.imi  tations 

This  study  is  restricted  to  classical  frequency  response  analysis. 
Therefore,  the  final  equations  are  linear  and  time  invariant.  Inputs 
to  the  system  are  restricted  to  lateral  stick  movements  and  turbulence. 
Since  the  F-15  aircraft  flight  controls  are  designed  to  maintain  coor- 
dinated flight  without  the  aid  oF  the  pilot,  rudder  inputs  are  not 
considered.  Also,  longitudinal  stick  inputs  as  a result  of  lateral 
stick  movement  have  not  been  considered  for  the  following  two  reasons: 

1.  The  aircraft  is  in  a trimmed  maneuver  requiring  no  back 
stick  pressure  to  maintain  the  desired  load  factor. 

2.  Precision  control  of  the  aircraft  load  factor  is  not  a 
requirement  for  this  system  (Ref  1:5.1). 

( The  effects  due  to  equipment  noise  or  measurement  errors  are  not  discussed 
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in  this  t'eport. 


In  order  to  simplify  the  equations  of  motion,  a constant  speed  of 
0.6  Mach  has  been  used  as  the  nominal  velocity.  For  a density  altitude 
of  p/p^  equal  to  0.971  (v/hicr.  has  also  been  assumed  to  be  a constant), 
the  total  steady  state  velocity  of  the  aircraft  is  665  feet  per  second. 
A constant  density  altitude  has  been  assumed  since  it  is  envisioned 
that  the  majority  of  the  weapon  delivery  modes  will  have  less  than  a 
15  degree  flight  path  angle  for  the  European  theater  where  low  overcast 
skies  will  restrict  the  height  of  the  aircraft  above  the  target.  For 
this  same  reason,  the  aircraft  trajectories  modelled  in  this  report  are 
either  level  or  shallow  descending  turns. 

I 1 


I 
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1 1 , Geometry  Equations 


The  linear  fire  control  bombing  law,  which  is  developed  in 
Chapter  3,  is  designed  to  null  the  perturbations  about  a nominal 
trajectory.  To  accomplish  this,  it  is  first  necessary  to  determine  the 
nominal  parameters  which  define  this  trajectory.  The  approach  used  in 
this  study  assumes  the  initial  conditions  are  the  nominal  values  which 
will  result  in  an  accurate  solution.  That  is,  if  no  disturbances  or 
inputs  are  added  to  the  system,  the  orientation  of  the  aircraft  is  such 
that  it  will  fly  to  an  accurate  release  point.  This  approach,  in  effect, 
freezes  the  aircraft  at  different  points  along  the  flight  path  so  that 
the  response  of  the  system  may  be  analyzed  at  these  various  conditions. 
The  equations  which  are  developed  in  this  chapter  determine  the  nominal 
conditions  for  the  angle  of  attack,  the  Euler  angles  and  the  different 
ranges  involved. 

First,  in  order  to  obtain  these  nominal  values,  it  has  been  assumed 
that  the  position  of  the  aircraft  relative  to  the  target  is  known.  This 
is  not  an  unreasonable  assumption  since  onboard  equipment  is  capable  of 
providing  this  information.  The  following  three  parameters  fix  the  air- 
craft position  relative  to  the  target: 

Rfo  - The  slant  range  from  the  aircraft  to  the  target 
h - Height  of  the  aircraft  above  the  target 
Xq  - Angle  between  the  velocity  vector  of  the  aircraft  and  the 
vector  from  the  aircraft  to  the  air  mass  aim  point,  a dis- 
tance Gq  above  the  target. 

Knowing  these  three  quantities  plus  the  velocity  of  the  aircraft  and 
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the  desired  time  of  fall  of  the  bomb,  the -necessary  nominal  values  for 
the  fire  control  law  can  be  determined.  These  nominal  values  are: 


1.  Gq  - Height  of  the  aim  point  above  the  target. 

2.  - Distance  from  the  release  point  to  the  aiming  point. 

3.  Y„  - The  angle  between  Rj  and  the  horizontal. 

4.  Rpjj  - Range  from  the  aircraft  to  the  aiming  point. 

5.  a - Angle  of  attack. 

O 

6.  iIiq.OqxJio  - Euler  angles  for  yaw,  pitch  and  roll. 

Although  the  nominal  time  of  fall  is  not  a fixed  parameter  in  the  actual 

"Firefly"  control  system,  it  has  been  given  a fixed  value  here  to 
simplify  the  calculations.  Where  an  infinite  number  of  solutions  existed 
before,  now  only  one  solution  exists. 

The  aimpoint  P is  defined  as  being  a gravity  drop  distance  Go 
directly  above  the  target  which  is  given  approximately  by: 


2 

Gq  ^ *5  Dp  9c 

'o 

(2-1) 

where  Dp  = 1 - Kn  (p/pq)  Vt  Tf 

0 0 

(2-2) 

Kjj  = Drag  constant  of  the  bomb 
6 1 

= 2.4  (10"  ) ft  ” for  low  drag  MK83 

p/pQ  = Relative  air  density 

= .971  for  this  analysis 

= Time  of  fall  of  the  weapon 
o 

V-r  = Inertial  velocity  of  the  aircraft 
’ 0 

= 665  ft/sec 
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Gravitational  constant 


9c  = 

= 32.2  ft/sec^  (Ref  2:14) 

The  release  range,  which  is  defined  as  the  distance  from  the  point 
of  release  of  the  bomb  to  the  aim  point  P,  is  given  by: 


Rr  = Op  (2-3) 

(Ref  2:15) 

Note  that  the  angle  between  the  horizontal  and  is  defined  as 

o 

positive  downward  as  shown  in  Fig.  3.  The  value  for  this  angle  is  given 
by: 


Y = sin"  (h/R,  ) (2-4) 

'o  I - 


The  angle  y is  also  defined  as  positive  downward  as  shown  in  Fig.  3. 
From  Fig.  3,  the  following  relationship  may  be  realized: 

R-r  sinT' 

Go=-I^ , (2-5) 

cos(Yo-Y  ) 

It  follows  that, 

GqCOS  Tq 

tan  Y (2-6) 

Rt  -G„sin  Y 
To  o 'o 


From  Fig.  3,  Rp  can  be  determined  using  the  law  of  cosines  and  the 

' O 

quadratic  equation: 


or. 


I 2 2 2 I L 

R-  = Ry  cosy  - (G  - Ry  sinT  y 

ro  In  O I- 


(2-7) 


(2-8) 
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where  Rp  is  the  position  vector  of  the  aircraft  relative  to  the  aim 

o 

point. 


Fig.  3.  Relationship  between  y'  and  Rp 

' 0 


The  next  unknown  which  must  be  determined  is  the  radius  of  turn. 

The  fire  control  bombing  law  is  based  upon  the  aircraft  flying  a circular 
flight  path  which  is  restricted  to  the  plane  containing  the  velocity 
vector  of  the  aircraft  and  the  vector  from  the  aircraft  to  the  air 

mass  aimpoint  (Rp^)-  The  vector  Rp^  is  determined  by  the  height  of  the 
aircraft  above  the  target,  the  time  of  fall  of  the  bonib,  which  in  turn 
determines  the  magnitude  of  G^,  and  the  range  from  the  aircraft  to  the 
target.  It  should  be  noted  that  the  vector  perpendicular  to  Rp^  which 
lies  in  this  plane  is  a horizontal  vector  and  is  designated  as  Rj^.  Since 
the  solution  is  based  upon  a constant  radius  turn,  the  load  factor  must 
be  adjusted  to  keep  the  aircraft  in  this  plane.  If  the  pilot  increases 
his  load  factor  without  increasing  his  bank  angle,  he  will  tilt  the  velo- 
city vector  of  the  aircraft  upward  out  of  the  present  plane.  This  will 
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result  in  a greater  time  of  fall  to  account  for  the  air  mass  aimpoint 

being  higher.  This  in  turn  raises  the  vector  L to  coincide  with  the 

'o 

new  plane  in  which  the  velocity  vector  now  rests.  In  this  way,  the 
actual  fire  control  system  is  able  to  update  itself  so  as  to  provide  a 
continuous  solution  as  the  aircraft  maneuvers. 


Fig.  4.  Relationship  Between  x and  R2 


From  Fig.  4,  and  using  the  law  of  cosines,  R,,  the  radius  of  turn  can 


be  defined  as: 


R,  = 


2 2 

- "h, 

2Rp^sfnx^ 


(2-9) 


Note  that  the  angle  is  defined  as  positive  for  a clockwise  direction 
and  is  negative  as  shown  in  Fig.  4.  Due  to  the  effect  of  the  sign  on 
different  angles,  all  equations  in  this  section  are  defined  for  an  air- 
craft attacking  the  target  from  an  orbit  which  is  counter  clockwise  with 
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Fig.  5.  Relationship  Between  Rp^  and  R^  to  the  x,  y,  z Coordinate  System 
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From  Fig.  5: 


Rp  cos  (y  -y') 
0 

Rp„  sin  iy-y') 

. r 0 0 


n = - sin  (y  -y')  i + cos  (y„~y')  k 


(2-11) 


Let  the  force  in  the  plane  which  produces  the  turn  be  defined  as: 


F = = m (A^  i + A j + k) 


Since  A and  n are  perpendicular: 


A’n  = - sin  (y  -y')  + A cos  (y  -y')  = 0 
A o ^0 

sin(Y  -y') 


cos(y^j-y') 


also:. 


A‘Rp  = A Rp  cos  (90°  -X  ) 

• rt  * rt 


= A^  Rp^  = A^  Rp^ 


A^  =*  A cos  (y  -y')  sin  X 

A 0 0 


Therefore: 


A^  = A sin  (yq^y' ) sin  x^ 


(2-12) 


(2-13) 

(2-14) 


(2-15) 

(2-16) 


(2-17) 


The  angle  between  the  horizontal  and  the  vector  A is  then  given  by: 


sinY 


* \ 


sin  Y*  = sin  (y  -y' ) sin  x 
0 o 


(2-18) 


* it 

It  should  be  noted  that  y is  positive  downward.  Knowing  y and 

Vt  ^ 

recalling  that  A can  be  set  equal  to — 1,  the  bank  angle  ()>,  can  now 

R2 

be  determined  with  the  aid  of  Fig.  6. 


\ it 

Fig.  6.  Bank  Angle  41  with  Respect  to  y 


attack  may  be  determined  as  follows: 


where 


“ *2  pVj  S Cj^ 


(2-22) 


pV-r 's  c, 

'o  a 


= .02974  radians  for  n = it 

For  any  given  flight  condition,  the  angle  of  attack  is  given  by: 


a 

o 


.02974  n 


.02974. 


To  • * 

9c  cosfo 


(2-23) 


In  order  to  determine  the  Euler  angle  g^,  it  is  necessary  to  find  the 
flight  path  angle  y"-  Let  the  velocity  vector  of  the  aircraft  be 
defined  as: 

i j + k (2-24) 

Since  the  velocity  vector  is  perpendicular  to  n.  the  dot  product  may 
be  written  as: 

V^«  n’  - sin  (y^-y')  + cos  (y^’y')  = 0 (2-25) 

tan  (y^-y')  (2-26) 


The  dot  product  for  the  velocity  vector  and  the  position  vector  yields: 


Vy.  Rp  “ Vt  Rp  cos  X 
*0  *^0  'o  *^0  0 

tFor  flight  conditions  given  in  Appendix  A. 


(2-27) 
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where  is  defined  in  Fig.'S. 


Vt^RPqCosx^  = U^Rp^  cos  (yq-y')  + W^Rp^  sin  (y^-y')  (2-28) 


Substituting  equation  (2-26)  for  W^,  the  longitudinal  velocity  is: 


Uo  = V.  cos  (y^-y')  cos 


(2-29) 


It  then  follows  that; 

W = Vt  sin  (y  -y')  cos  X (2-30) 

O I o O o 

The  flight  path  angle  can  now  be  determined  by: 

W 

sin  y"  = — = sin  (y  -y')  cos  x (2-31) 

woo 

''To 

also.  V/  = U/+  V/+  W/  (2-32) 

o 

''o  " [''T^  V-r"cos''(Y^-Y')cosV  VT%in"(Yo-Y')cos"xJ'' 

= (1-cos^Xq)^ 


Vt  sin  X 

To  o 


(2-33) 


To  describe  vectors  in  the  aircraft  reference  frame,  it  is 
helpful  to  use  the  Euler  transformation  matrices  as  follows: 


X 

y 

A 

Z 


(2-34) 


Where  x,  y and  z describe  the  body  fixed  aircraft  coordinate  system 

A A A 

and  i,  j and  k relate  to  the  earth  fixed  coordinate  system. 
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I 


1 0 0 

0 cosiji^  sin>|i^ 

0 -siriij  cosii 
o o 


coso^ 

0 

-sine 

0 

1 

0 

sine 

0 

0 

cose 

1 

cosily 

sin^,^ 

0 

-sin* 

0 

COSi}.^ 

0 

0 

0 

1 

cosii  cose 

o 

cosi>^sin0^sini{)^-sinip^cos.{)^ 

cosij/^sinOQCOSij)^+sinij;^sin(i)^ 


sinii)  cose 
o o 

s 1 nii/^s  1 ne^s  1 p<ji^+cos(|;^cos(()^ 
s i nij;^s  i n 0^  cos  - cos  s i riifiQ 


The  aircraft  velocity  may  then  be  expressed  as: 


Vo 

= 

V-r  COSa„ 
‘o  ° 

0 

■w 

Vj  COs(y  -y') 
'o  o 

0 

W 

0 

V-r  sinao 
. 0 

Vt  Sin(yQ-y') 

. 0 J 

where  V-  is  aligned  with  Rr 


-sine 

0 

cose  sinA 
0 

COSO  COSd) 

O 0 


(2-35) 


For  level  flight,  y"  = Y - y'  = 0»  the  above  velocity  equations  then 

o 

reduce  to: 


coso^  = cosiji^cose^  (2-36) 

0 = cosi^/^sine^sini)!^  - sinij»^coS({i^  (2-37) 

sina  = cosii)^sine^coS(j>^  + sin4/^sin(^^  (2-38) 

The  angles  and  e^  are  unknown  in  the  above  equations.  The  functions 

containing  may  be  replaced  by  functions  containing  a^,  e^,  and 

It  is  then  possible  to  solve  for  o as: 

O 


( 

I 


tane  = tana  cos<4 
0 0 O 

For  a climbing  or  descending  turn, 

tan(0Q+Y")  = tanaQC0S(>^ 


(2-39) 


(2-40) 
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tanOg  = 


tann^  cos<(>^^  - tan  y" 

1 + tann  cosiji  tan  y” 
0 o 


(2-41) 


The  last  remaining  unknown  is  the  Euler  angle  Using  the 

Euler  transformation  matrix  and  the  more  general  case  where  Vt 

* o 

is  not  aligned  with  Rn  , the  aircraft  velocities  may  be  expressed  as: 

» Q 


Vt  COSa 
‘ 0 


Vt  sina„ 
. To 


Mj  COS(yg-y')cOSX^ 

[t] 

V-p  sin(y^-Y')cosx^ 

- * o 


(2-42) 


Where  the  velocities  for  the  earth  fixed  coordinate  system  cone  from 
equations  (2-29),  (2-30)  and  (2-33).  The  only  unknown  in  this 
expression  is  ip^which  may  be  written  as: 

COS4<q=  ^COSaQCOs(y^j-y  ' )cOSX^COS0^  + COS^i^sina^COS  (y^-y  ' )cOSXQSi  n0^  - 

siniji^sina^sinx^  J / [cos^(y^-y ' )cos  Xq  + sin  xj  (2-43) 
Table  I presents  a summary  of  the  geometry  equations  defining  the 


steady  state  input  parameters  used  in  the  fire  control  equation. 
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Table  I 


Summary  of  Geometry  Equations 

Fixed 


III.  Fire  Control  Law 


As  mentioned  in  the  introduction,  the  fire  control  bombing  law 
controls  the  bank  angle  of  the  aircraft  in  order  to  achieve  a solution. 
It  is  able  to  accomplish  this  by  maintaining  a constant  load  factor 
while  varying  the  bank  angle  v/hich  raises  or  lowers  the  velocity 
vector.  When  the  velocity  vector  passes  above  the  target,  it  should  be 
at  the  correct  inclination  so  that  the  bomb  will  cover  the  horizontal 
distance  to  the  target  during  the  time  it  takes  the  bomb  to  fall  to 
the  target.  For  the  purposes  of  this  study,  the  fire  control  law  can 
be  broken  dov/n  into  tv/o  segments,  the  steady  state  condition  portion 
which  represents  the  nominal  solution  and  the  perturbed  condition 
v/hich  represents  the  errors  due  to  external  disturbances. 
fionl inear  Fire  Control  Law 

It  is  first  necessary  to  express  the  fire  control  bombing  law  in 
terms  of  the  geometry  parameters  developed  in  chapter  two.  The  nominal 
flight  path  of  the  aircraft  has  been  determined  by  fixing;  the  position 
of  the  aircraft  relative  to  the  target,  the  height  of  the  air  mass 
aimpoint  above  the  target,  and  the  velocity  of  the  aircraft.  With 
the  geometry  thus  defined,  it  is  now  possible  to  write  the  nonlinear 
fire  control  law  in  terms  of: 

1.  The  state  variables:  a,  v,  0,  <». 

2.  The  position  vectors:  G,  Rp,  Rj^,  Rj,  h 

3.  The  parameters:  T^,  y*  Vj  , I 
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where 


(3-1) 


G.  .5[l  - K„  (p/„^)  Vj  Tf]  t/ 

Rr  = [l  - Kg  (p/»o>  "t  h]  ''t  T, 


_l/  G cos  Y \ 

cos 

tan  ( ) 

- 

yRl--G  sinyj 

I 

. 2 
Sin 


tan 


_i/  G cos  Y 
^R^-G  siny/J 


(3-2) 

(3-3) 

(3-4) 


Since  G,  Rj^,  Rp  and  y can  be  expressed  in  terms  of  V-j-,  T^,  h and  Rj, 
the  variables  of  the  system  are  then:  a,  S',  o,  4>,  Vj,  T|:,  h,  Rj  and  x. 
Note  that  V-j-,  T^,  h,  Rj  and  x are  the  parameters  which  have  been 
fixed  in  order  to  set  the  trajectory  of  the  aircraft. 

Using  the  notation  in  the  "Firefly"  manual,  the  fire  control 
bombing  equation  is  written  as: 


P 


c 


(3-5) 

. (Ref  2:26) 


where 


= Roll  rate  command 


K = System  gain 

The  remaining  terms  may  be  written  in  terms  of  the  gectnetry  inputs 

previously  mentioned.  Cp  is  defined  in  reference  2 as: 

w 


, , (3-6) 

(Ref  2:21) 


where 


S = coso 
''u 


(Ref  2:20) 


(3-7) 
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V V 


(3-8) 


(Ref  2:14) 

Note  that  u,  v and  w correspond  to  the  x,  y,  and  2 body  axes  of  the  air- 
craft. In  order  to  simplify  the  above  expressions,  the  ejection  velcrcity 
of  the  bomb,  the  target  velocity,  and  the  steady  wind  velocity  have  all 
been  assumed  to  be  zero.  With  the  aid  of  the  Euler  transformation  matrix 
as  defined  in  chapter  two,  R-j-  and  G may  be  written  as: 

= R.j.|^cosy  (cos'i'sinosin'i’  - sin'i'cos'i')  + sinycosOsin'fj  (3-9) 


Gy  = G cososini' 


(3-10) 


is  expressed  as: 


2 2 

Rp  - Rp 


(3-11) 
(Ref  2:26) 


and  needs  no  further  defining  since  it  is  already  in  terms  of  the 
geometry  inputs. 


The  term  n is  defined  as: 
w 


^ (3-12) 

Vj 

' (Ref  2:22) 

v/here  A^,  the  acceleration  in  the  lateral  direction  may  be  written  as: 


'"a  ''t 

A = — 1 + g + - 

m ^v  m 


(3-13) 


Since  the  engines  are  aligned  along  the  x axis, 


(3-14) 
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For  the  steady  state  flight  condition  where  the  aircraft  is  assumed 
to  be  in  coordinated  flight  with  zero  sideslip: 


= g cose  sin(j) 

C 0 0 


(3-15) 


The  last  undefined  term  in  the  roll  rate  command  equation  if 
which  is  expressed  as: 


= Rp  sin  X (3-16) 

(Ref  2:19) 

It  should  be  noted  that  this  fire  control  bombing  law  contains  the 
error  for  the  nominal  conditions.  This  error  may  be  expressed  as; 

Error  = (3-17) 

The  terms  defining  the  nonlinear  fire  control  equation  are 
summarized  in  Table  II. 

Linear  Fire  Control  Law 

In  order  to  study  the  effects  of  the  perturbations  on  the  system, 
the  equations  can  be  linearized  so  that  frequency  response  analysis 
may  be  used.  To  accomplish  this,  the  variables  which  are  responsible 
for  causing  a change  in  the  roll  rate  command  must  be  determined.  The 
perturbations  of  interest  are  the  ones  which  change  the  velocity  vec- 
tor. These  are:  u,  0,  a,  p,  r,  i|),  o,  41,  (S^y  and  6j^.  Also,  as  these 
state  variables  change,  they  produce  a change  in  the  time  of  fall  and 
the  heading  angle  x.  These  changes  may  be  denoted  as  ty  and  x respect 
ively.  Since  Ry  is  of  a large  magnitude,  even  close  to  the  release 
point,  it  has  been  assumed  to  be  a constant.  Also,  due  to  the  large 


25 


values  of  and  h,  y has  also  been  assumed  to  be  unaffected  by  per- 
turbations. Finally,  the  magnitude  of  Rp  is  independent  of  the  time 
of  fall  to  a first  order  of  approximation.  (Ref  2:17). 

As  and  x cannot  be  considered  constants,  it  is  necessary  to 
define  these  parameters  since  their  perturbations  will  have  an  effect 
on  the  solution.  The  equation  which  governs  the  time  of  fall  of  the 
bomb  for  a vacuum  is: 


h = Z + ->5  g^  t/ 


(3-18) 


where  Z is  the  initial  vertical  velocity  in  the  earth  fixed  coordinate 
reference  frame.  Using  the  Taylor’s  series  expansion  and  retaining 
only  the  linear  terms  yields: 


Ah  = z T^  + Z^  tf  + Tf^  tf 


Since  we  assume  Ah  is  equal  to  zero,  t^^  may  be  written  as: 


- z Tf 

t.  = . 


v/here  Z is  determined  from  the  following  relationship: 


“ ■ 

- - 

X 

u 

Y 

= N" 

Vt  sin  i 
•o 

Z^ 

Vt  sin  a 

L *0  J 

or,  Z = - U sine  + Vt  sino„  cose  cos<p 

O O O In  o 0 


(3-19) 


(3-20) 


(3-21) 


(3-22) 


Again,  using  the  Taylor's  series  expansion  the  linearized  equation  for 
z is : 
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(3-23) 


z = u (-sino^)  + b(Vi-  coso^siniji^)  + u(W^  cose^cosiji^)  - 
o{U  cose  + W sine  cosi}>  ) - ifi(W  cose  sina.  ) 

0 O 0 0 *^0  ' 0 o 

The  angle  A is  defined  as  the  angle  between  the  velocity  vector  and 
Rp,  or: 

• Rp  = Rp  cos  A (3-24) 

The  roll  rate  command  equation  may  now  be  linearized  by  taking  the 
Taylor's  series  expansion  of  the  equations  for:  T,,  A,  Cp  , Oi-[^, 

Rq^.  and  Pj,.  Their  respective  values  for  each  of  the  perturbed  varia- 
bles are  defined  in  Fig.  21.  (Elements  (14,1)  through  (20,5)). 

The  roll  rate  command  is  then  passed  through  an  integrator  to 
obtain  a bank  angle  steering  bar.  This  step  is  accomplished  since  the 
literature  refers  to  roll  task  models  where  the  bank  angle  is  being 
controlled  and  not  the  roll  rate. 
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Table  II  (Continued) 


Summary  of  Nonlinear  Fire  Control  Law 


Fixed 

Parameters 

jt 

:e 

Variables 


I V , Turbulence  Rofiponse  Model 


To  determine  the  aircraft  response  to  wind  gusts,  it  is  necessary 
to  represent  the  atmospheric  turbulence  by  nathernatical  models.  This 
can  be  done  by  representing  the  wind  gusts  as  input  power  spectral 
densities  (PSD).  It  is  assumed  here  that  the  turbulence  is  isotropic, 
homogeneous  and  of  a stationary  nature  for  altitudes  above  1750  feet 
AGL.  (Ref  3:419-421).  The  input  PSD  can  then  be  represented  as: 

PSOi^  = <^i  = |T.  (co)I^D  (4-1) 

g g 

where  T^.  (w)  is  the  transfer  function  describing  the  shaping  filter  which 

g 

transforms  white  noise  of  intensity  D into  colored  noise  representing  the 
wind  gusts  for  that  particular  direction.  The  PSDs  for  the  wind  gusts 
along  the  lateral  and  vertical  axes  are  represented  by  (j>  and  while 

g g 

the  PSDs  for  the  wind  gusts  angular  velocities  about  the  roll,  pitch,  and 
yaw  axes  are  represented  by  Ppp 

analysis,  the  intensity  or  rms  value  of  the  white  noise  (D)  may  be  assumed 
to  be  unity. 

The  Oryden  spectral  form  of  the  turbulence  will  be  used  in  this 
analysis  since  it  is  a simpler  model.  For  flying  qualities  analysis, 
the  Dryden  model  yields  essentially  the  same  results  as  the  more 
complicated  VonKarman  model  (Ref  3:422).  The  Dryden  form  for  the 
continuous  random  nwdels  are: 
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2 / "‘-W  \ 1/3 
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Vt  L„  iWib. 
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(4-4) 


to  <{ly;  ( to ) 

t’q  ^ “Tf — h 

Vt  r 1 +/4b 


‘ 1^)  1 


(4-5) 


(ir  (w) 

g 


to  tpv  (to) 

Vt^  r 1 +/3b  0) 

Uv 

L 1 0 


(4-6) 

(Ref  3:459,460) 


v/here 


= The  root-mean  square  intensity  of  the  wind  gusts. 


L.  = Scaling  factor  for  the  wind  gust. 


to  = Frequency  in  radians/second. 


b = Wing  span  of  the  aircraft. 


Since  the  has  little  effect  on  high  speed  aircraft  due  to 
their  high  momentum,  is  not  considered  in  this  analysis. 

For  the  clear  air  turbulence  nwdel , the  scales  for  the  Dryden 


form  are; 


Ly  = Ly  = = 1750  feet 


(4-7) 
(Ref  3:444) 


for  altitudes  above  1750  feet  above  the  terrain.  The  root-mean-square 


intensity  is  defined  as  approximately  six  feet  per  second  for 
altitudes  around  2500  feet  above  the  terrain.  The  intensities 
and  ow  may  then  be  obtained  from  the  following  relationship: 


2 2 2 
*-U  *"7 


(4-8) 
(Ref  3:435) 


This  results  in  the  three  intensities  being  equal  or  = 6 feet/sccond. 
Noting  that  and  Wg  are  correlated,  it  is  possible  to  represent 
^ as: 


^q  " l'''qjjw)  I ‘ V 
^g  y 9 


(4-9) 


where 


(s)  = 


1 -H  4bs 
ttVt 


(4-10) 


Likewise,  <>  and  (}>  are  correlated,  where  (fr  represented  as: 


1>^  (u)  = IT^  (ju)  I • 4>y  (w) 


(4-11) 


where 


T^  (s)  = 
9 


1 + 3bs 
ttVt 

* m 


(4-12) 
(Ref  3:459-460) 


Since  the  equations  describing  the  basic  aircraft  are  in  terms  of  a 
C..d  Bi  it  should  be  noted  that: 


♦.g  <“>  = i'agl  ' 


(4-13) 


> 


(4-14) 


•!>  (w) 


where 


(4-15) 
(Ref  3:421) 


The  three  power  spectral  densities,  » snd  <p  , which  are  the 

'^g  '^g  ^g 

uncorrelated  inputs  to  the  system  are  plotted  in  figures  7 and  8 . 
Since  is  equal  to  and  is  equal  to  L^,  is  equal  to  . 

To  analyze  the  response  of  a system  to  an  input  PSD,  it  is  necess- 
ary to  obtain  the  output  PSD  which  is: 


PSDq^^  = |G  (ja,)|"  • PSD.^  (4-16) 

where  G (ju)  is  the  open  loop  transfer  function  of  the  system  being 
analyzed.  G (jw)  is  determined  by  replacing  all  the  linear  and 
angular  velocities  in  the  aerodynamic  terms  of  the  aircraft  equations 
with  the  following: 


Ua  = u + Ug  Pa  = P ^ Pg 

0A  = e + Sg  qA  = q + qg 

aA  = a + Og  '"a  ""g 

These  wind  gust  components  are  then  treated  as  inputs  to  the  system. 

In  matrix  form,  the  system  equations  may  be  written  as: 

[Bjgg*[B5]rg  (4-17) 

where  ^A  j - Matrix  which  defines  the  system  components. 
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- The  vector  representing  perturbation  variables. 

- The  vector  which  represents  the  inputs  to  the  system. 
(See  Appendix  D for  a more  complete  explanation). 

For  a wind  gust  Bol  is  zero.  Also,  since  and  r„,  and  a and 

L H g 9 g 

are  correlated,  they  may  be  combined  as: 


[b]=[b, 


\ * t 


Be 


+ Bl 


J 


(4-18) 


for  the  power  spectral  density  inputs  <ti  , <{>  , and  <j)  respectively.  To 

g g *^9 

simplify  the  description  for  these  combined  inputs,  the  wind  gusts  which 
are  related  to  the  <f>  inputs  will  be  called  lateral  gusts,  while  the 

g g 

inputs  will  be  referred  to  as  longitudinal  gusts  and  the  ij>  inputs  will 
be  annotated  as  P gusts  or  roll  gusts.  The  combined  effect  of  these  three 
gusts  will  be  referred  to  as  a composite  gust.  The  value  of  the  composite 
gust  may  be  calculated  by  summing  the  squared  values  for  each  of  the  three 
individual  gusts  and  then  taking  the  square  root  of  the  sum.  Figure  23  in 
Appendix  D presents  the  B matrix  as  defined  in  the  computer  program  used  for 
this  study.  These  wind  gust  inputs  originate  in  the  force  and  moment 
equations  of  the  aircraft  and  the  flight  control  equations  where  f. 

y 

and  f.  are  fed  back  into  the  pitch  and  yaw  axis  equations.  Figure  24  in 
''z 

Appendix  D represents  the  B matrix  where  the  flight  controls  have  been 


modified  for  the  automatic  fire  control  system.  The  P parameter  which 
controls  the  sequence  of  the  inputs  is  also  defined  in  Appendix  D. 


PSD 


V.  Pilot  Model 


:1 


Background 

In  recent  years  a groat  deal  of  work  has  been  done  in  an  attempt 
to  represent  the  human  pilot  by  an  analytical  model.  These  pilot 
models  are  used  in  conjunction  with  system  models  to  form  predictions 
or  explain  the  behavior  of  the  closed  loop  system  where  the  pilot 
closes  the  loop.  Systems  Technology,  Inc.  has  been  prominent  in 
developing  different  pilot  models  and  the  reader  is  referred  to 
reference  4 for  a thorough  discussion  of  the  topic.  The  general  model 
developed  in  reference  4 is  given  in  its  simplified  form  as: 


where 


jto  + 1 
ju)  + 1 


T,,  j(i)  + 1 


(5-1) 
(Ref  4:17) 


Kp  = Pilot  gain 


g-jioT  _ jpapsport  delay  due  to  basic  latencies  such  as 
nerve  condition,  data  processing,  etc. 


\Tj  jo,  + 1/ 


= Pilot  adjustment  parameters  and  equalization 
characteristics  the  pilot  will  adopt  to  suit 
the  task. 

= First  order  neuromuscular  lag  terms. 


This  model  can  be  simplified  still  further  by  combining  the  transport 
delay, due  to  basic  latencies  with  the  neuromuscular  system  lag  term. 
This  combination  is  often  referred  to  as  the  effective  time  delay,  Xg. 
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These  simplifications  lead  to  the  more  commonly  used  model: 


(5-2) 
(Ref  5:234) 


where  ju  = s for  continuous  random  like  inputs. 

Crossover  Model 

The  pilot-vehicle  models  are  most  useful  for  compensatory  track- 
ing tasks  where  the  pilot  acts  upon  a displayed  error  due  to  some  random 
input.  The  transfer  function  describing  the  controlled  element  is 
Yj.  (s)  while  the  pilot's  control  action  which  is  linearly  correlated 
with  the  input  is  described  by  Yp  (s).  In  order  to  minimize  the  error, 
the  amplitude  ratio  of  the  open  loop  frequency  response,  IYolI  = Vc  • 
should  be  very  large  over  the  frequency  of  the  input  band  width  and  very 
small  outside  this  range.  (Ref  5:233). 

The  pilot  attempts  to  achieve  this  condition  by  adopting  suffici- 
ent lead  or  lag  equalization  so  that  the  slope  of  |Yqj^|  lies  very  close 
to  - 20  dB/decade  in  the  region  of  crossover  frequency,  This  crossover 
frequency  occurs  where  |Yqj^|  = 1.0  or  0 dB.  For  small  tracking  errors, 
should  be  greater  than  the  input  frequency,  . The  pilot  can  adjust 
this  crossover  frequency  by  adjusting  his  gain,  Kp.  This  leads  to  the 
"two  parameter  crossover  model  which  can  account  for  most  of  the  signifi- 
cant open  loop  data  trends  in  the  important  crossover  frequency  region. 

The  crossover  model  is: 

u> 

% <*'  " V ''c  ■ 

s 
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near  (Ref  5:234).  Here,  is  equivalent  to  the  pilot  or  loop 

gain  while  represents  the  lead  or  lag  induced  by  the  pilot.  If  the 
form  of  is  knov/n,  Yp  can  be  determined  from: 


-T  s 
_a)ce  e 


H s Y^ 

In  reference  5,  typical  aircraft  control  tasks  have  been  determined 
where: 


(5-4) 


Yc=- 


s(Ts  + 1) 

represents  the  roll  angle  being  controlled  by  lateral  stick  inputs 

(Ref  5:236).  Y can  then  be  written  as: 

P 

Yp  = Kp  (Tl  s + 1)  e-V 


(5-5) 


(5-6) 


This  indicates  that  a typical  system  would  show  a - 40  dB/decade  slope 
in  the  region  prior  to  the  crossover  frequency.  The  pilot  then  adds 
enough  lead  to  make  the  slope  in  this  region  - 20  dB/decade. 

In  reference  6,  an  analysis  was  conducted  where  the  roll  task  was 
studied  with  regards  to  pilot  models.  In  this  same  report,  it  was  deter- 
mined that  the  pilot  model  could  be  represented  as: 


Yp  - Kp  (.5  s ♦ l)e' 


'.3$ 


(5-7) 


for  the  roll  task  which  consisted  of  attempting  to  maintain  a fixed  bank 
angle  in  a turbulent  field  (Ref  6.12). 

The  pilot  model  represented  in  equation  (5-7)  is  the  model  used  in  this 
thesis.  The  reasons  for  this  choice  are: 

1.  The  same  task  is  being  accomplished. 

2.  The  Inputs  to  the  system  are  modeled  as  turbulence. 
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The  pilot  gain  is  determined  by  placing  the  crossover  frequency  at 
approximately  four  radians  per  second  for  an  open  loop  controlled 
system  with  a - 40  dB  slope  prior  to  the  crossover  point  (Ref  5:237). 
System  Application 

The  pilot  model  can  be  added  to  the  system  as  shown  in  figure  9. 


Fig.  9.  Closed  Loop  System  with  Pilot  Model 


The  reference  angle  (|)^  may  be  set  equal  to  zero  since  41^  represents  the 
perturbations  about  this  reference  bank  angle.  With  this  simplification 
the  pilot's  horizontal  stick  inputs  can  be  written  as: 


'"as  " 


Kp  {.5s  + l)e‘ 


:c 

s 


(5-8) 


where  is  equal  to  p^/s. 
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With  these  inputs,  the  primary  roll  axis  equations  in  Appendix  C can 
now  be  written  as: 


6^  (23.87)  + 


= 0 


(5-9) 


p p • 3 , /.0019s  + .0074\ 

p(.0625)  + 6p.^(2)  + ^ ['<p(-5s  + l)e  ' I = 0 (5-10) 


For  the  basic  aircraft,  the  pilot  model  will  be  attempting  to  null  the 
perturbation  roll  angle  <j).  The  primary  roll  axis  equations  will  be  the 
same  as  the  above  except  that^  is  replaced  by  <j). 
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VI.  System  Evaluation  and  Results 

The  preceeding  chapters  along  with  the  appendices  develop  the 
equations  which  define  the  system  and  describe  the  computer  program  to 
solve  for  the  frequency  response.  The  purpose  of  this  chapter  is  to 
evaluate  those  areas  which  will  indicate  how  the  system  responds  to 
random  inputs  and  to  compare  these  responses  for  different  configurations 
of  the  system. 

Form  of 

The  first  item  to  be  considered  is  the  form  of  the  transfer 
function  for  the  complete  system  including  the  fire  control  system. 

Figure  10  shows  the  Bode  plot  asymptotes  for  the  basic  aircraft  and  the 
complete  system  where  the  input  is  an  aileron  deflection  of  unit  magnitude. 

The  output  is  41^  for  the  complete  system  and  <t>  for  the  basic  aircraft. 

Recalling  from  chapter  5,  the  pilot  must  adjust  his  lead  or  lag  to 
achieve  a -20  dB/decade  slope  for  the  combined  Yp  Y^  response  for  a good 
portion  of  the  region  prior  to  the  crossover  frequency.  He  then  adjusts 
his  gain  to  place  the  crossover  frequency  somewhere  around  4 radians/ 
second.  For  the  basic  aircraft,  all  the  pilot  needs  to  do  is  adjust  his 
gain  since  this  system,  which  already  demonstrates  a -20  dB  slope,  needs 
no  further  equalization.  For  the  complete  system  the  pilot  must  not 
only  adjust  his  gain  but  he  must  also  add  some  lead  to  reduce  the 
-40  dB/decade  slope  to  approximately  -20  dB/decade.  The  pilot  model 
which  was  chosen  in  chapter  5 is  capable  of  producing  the  desired 

f 

slope  of  -20  dB/decade  by  adding  the  .5  seconds  of  lead  time.  This  only  i 

f 
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leaves  the  pilot  gain  to  be  adjusted  so  as  to  achieve  the  desired  cross 
over  frequency. 


Fig.  10.  System  Response  to  a Unit  Lateral  Input 


Pilot  Gain 

The  pilot  gain  for  this  system  is  measured  as  pounds/radian. 

As  the  pilot  gain  is  increased,  the  crossover  frequency  increases 
resulting  in  a larger  bandwidth.  The  pilot  gain  for  the  different 
configurations  and  different  trajectories  is  determined  by  increasing 
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the  open  loop  system  gain  to  achieve  a crossover  frequency  of  approxi 
mately  4 radians/sccond.  The  pilot  then  adopts  this  gain  so  as  to 
maintain  the  sanK  crossover  frequency  for  the  closed  loop  system. 
Table  III  shows  the  required  pilot  gains  to  achieve  a crossover 
frequency  of  approximately  4 radians/second.  It  should  be  noted  that 
for  the  basic  aircraft,  the  pilot  gain  is  independent  of  bank  angle. 
However,  for  the  case  where  the  pilot  is  following  the  command  bank 
angles  as  displayed  on  the  HUD,  the  pilot  gain  increases  as  the  bank 
angle  increases.  This  characteristic  does  not  have  any  significant 
influence  on  the  pilot  rating,  since  it  is  easily  compensated  for  by 
the  pilot. 


Table  III 
Pilot  Gain 


Flight  Condition 
Straight  and  level 
Level  turn 
Descending  turn 


for  (j> 


for  ({ic 


25  5 

25  80 

25  80 

♦ (Uj.  = 4 rad/sec 


RMS  Values  for  and  ^ 

In  order  to  analyze  the  bank  angle  system  response,  three  types 
of  trajectories  were  considered: 

1.  Straight  and  level  flight  (the  aircraft  heading  remains 
aligned  with  the  target  while  maintaining  a constant  altitude 
of  2500  feet  above  the  target). 
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2.  Level  turn  (the  initial  aircraft  heading  is  30°  off  the 
target  at  a slant  range  of  12,000  feet  while  maintaining  a 
constant  altitude  of  2500  feet  above  the  torget). 

3.  Descending  turn  (the  initial  aircraft  heading  is  again  30° 
off  the  target  at  a slant  range  of  12,000  feet  and  in  a shallow 
descent  of  approximately  10°). 

With  random  wind  gusts  as  inputs,  each  of  these  trajectories  are  evalua- 
ted to  determine  the  rms  values  for  and  This  is  done  for  the 
following  configurations : 

1.  Basic  aircraft  (no  pilot  and  no  fire  control  system),  see 
Table  IVa. 

2.  Basic  aircraft  with  pilot  closing  the  loop  by  trying  to  null 
the  i{)  perturbations,  see  Table  IVb. 

3.  Aircraft  plus  fire  control  system  (no  pilot),  see  Table  Va. 

4.  Aircraft  plus  fire  control  system  with  pilot  closing  the  loop 
by  trying  to  null  the  41^,  perturbations,  see  Table  Vb  and  Vc. 

5.  Fire  control  system  output  signal  fed  directly  to  the  F-15 
flight  controls,  see  Table  Via  and  VIb. 

Tables  IV,  V and  VI  summarize  the  rms  values  for  <t>  and  41^  which  are  in 
degrees.  Comparing  Table  IVaand  Table  Vc,it  is  interesting  to  note  that 
in  two  of  the  three  cases,  the  perturbed  bank  angle  rms  values  for  the 
closed  loop  system  actually  increased  over  the  open  loop  system  even 
though  the  pilot  v;as  able  to  essentially  null  the  perturbed  error 
signal.  Table  IV  indicates  that  the  pilot  is  much  more  successful  at 
reducing  the  rms  perturbed  bank  angle  values  by  trying  to  null  4 instead 
of  41c.  When  for  the  automatic  system  is  compared  to  for  the 

manual  system,  it  is  obvious  that  the  pilot  is  better  able  to  null  the 


) 
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error  signal  than  the  flight  controls.  However,  the  perturbed  bank 
angle  is  generally  greater  for  the  pilot  than  the  autopilot.  As  a 
final  observation,  it  can  be  noted  that  the  open  loop  ox  is  extremely 
large  inferring  the  system  will  direct  large  bank  angle  commands  if  'the 
perturbations  are  allowed  to  go  unattended.  The  probable  cause  for  the 
high  degree  of  sensitivity  for.  the  straight  and  level  case  is  the 
denominator  of  the  bombing  equation  (Eq.  3-5)  becomes  very  small  as  X 
goes  to  zero. 

Table  IV 

Bank  Angle  Response  for  Aircraft  without  Fire  Control  System 


a.  a*  " RMS  Value  for  Open  Loop  Perturbed  Bank  Angle 
Type  Gust  Straight  & Level  Level  Turn  Descending  Turn 


I Lateral 

1.33 

1.70 

1.47 

j Vertical 

.01 

.21 

.13 

P 

1.33 

1.67 

1.63 

; Composite 

1.93 

2.41 

2.27 

b. 

- RMS  Value  for  Closed  Loop  Perturbed  Bank  Angle 

1 

i Type  Gust 

Straight  & Level 

Level  Turn 

Descending  Turn 

Lateral 

.43 

.57 

.49 

' Vertical 

1 

.0 

.01 

.01 

' P 

1 

.18 

.19 

.19 

1 Composite 

.47 

.61 

.53 

RMS  Values  are  in  Degrees. 


46 


Table  V 


Bank  Angle  Response  for  Aircraft 

with  Fire  Control  System  and  Pilot 

3.  Oa 

•fcoL 

RMS  Value  for  Open  Loop  Perturbed  Command  Bank  Angle 

Type  Gust 

Straight  & Level 

Level  Turn 

Descending  Turn 

Lateral 

Li  mi  t 

7.96 

14.16 

Vertical 

3.44 

3.84 

3.33 

P 

Limit 

34.44 

30.25 

Composite 

Limit 

39.03 

38.86 

‘^CCL 
Type  Gust 

RMS  Value  for  Closed 

Straight  & Level 

Loop  Perturbed 

Level  Turn 

Command  Bank  Angle 

Descending  Turn 

Lateral 

5.12 

.03 

.11 

Vertical 

.01 

.01 

.01 

P 

.64 

.04 

.04 

Composite, 

5.15 

.05 

.11 

c.  0,  - RMS  Value  for  Closed 

■^CL 

Type  Gust  Straight  & Level 

Loop  Perturbed 

Level  Turn 

Bank  Angle 

Descending  Turn 

Lateral 

5.83 

2.19 

3.16 

Vertical 

.01 

.15 

.11 

P 

.08 

.09 

.09 

Composite 

5.84 

2.19 

3.16 

RMS  Values  are  in  Degrees. 
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Table  VI 

Bank  Angle  Response  where  Fire  Control  System 
is  Connected  to  the  Flight  Controls 

a.  ax  - RMS  Value  for  Perturbed  Command  Bank  Angle 
Type  Gust  Straight  & Level  Level  Turn  Descending  Turn 


Lateral 

6.02 

2.89 

3.28 

Vertical 

.01 

.19 

.19 

P 

.3 

.25 

.24 

Composite 

6.05 

2.95 

3.36 

b. 

- RMS  Value  for 

Perturbed  Bank 

Angle 

Type  Gust 

Straight  & Level 

Level  Turn 

Descending 

Lateral 

4.81 

2.17 

2.76 

Vertical 

.01 

.17 

.13 

P 

.10 

.26 

.25 

j Composite 

4.83 

2.19 

2.78 

Error  Analysis 

The  distance  the  bomb  misses  the  target  is  an  excellent  tneasure 
for  determining  if  a bombing  system  is  of  any  use.  Since  the  solution 
for  this  fire  control  system  is  being  continually  updated,  only  the 
perturbations  injected  into  the  system  at  the  time  of  release  will 
result  in  an  error.  This  error  or  miss  distance  is  related  to  the 
difference  between  the  magnitude  and  orientation  of  the  perturbed 
velocity  vector  and  the  nominal  velocity  vector.  Since  the  miss 
distance  is  measured  in  the  earth  fixed  coordinate  system,  it  is 
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advantageous  to  use  the  Euler  transformation  matrix  as  follows: 


X 

u 

. 

r 1 

Y 

Vt  sin3 
' 0 

Z 

V-r  sina 

L 

L 'o  J 

where  X,  Y and  Z are  total  velocities  in  the  earth  fixed  reference  frame. 
The  X,  Y and  Z equations  may  be  written  as: 

X = U (cos'i'  cos©)  + Vj  sine  (cosf  sin©  sin4i  - sin*}*  cos4)  + 

0 

Vj  sina  (cosy  sin©  cos4>  + siny  sin4>)  (6-1) 

0 

Y = U (siny  cos©)  + Vj  sine  (siny  sin©  sin4>  i-  cosy  cos4>  ) + 

Vj  sina  (siny  sin©  cosit  - cosy  sin4<)  (6-2) 

Z = U (-sin©)  + Vj  sine  (cos©  sin<i))  + Vj  sina  (cos©  cos'i')  (6-3) 

o o 

where  the  above  variables  include  the  steady  state  and  perturbed  quan- 
tities. By  using  the  Taylor's  series  expansion  and  retaining  only  the 
linear  terms,  the  expressions  defining  the  perturbation  values  of  x,  y 
and  z may  be  determined. 

The  lateral  error  may  then  be  written  as 

e = y (t^  + t^)  (6-4) 

and  the  longitudinal  error  may  be  written  as 

Neglecting  the  products  of  the  perturbations,  the  lateral  and  longitud- 
inal errors  may  be  written  as; 


(6-7) 


The  total  error  may  be  written  as: 

er  = (e/  (6-9) 

Again  using  the  random  wind  gusts  as  inputs,  power  spectral  density 
analysis  can  be  performed  to  calculate  a rms  miss  distance  based  upon 
the  preceeding  equations.  Tables  VII,  VIII  and  IX  show  the  average  miss 
distances  for  the  diffet'ent  trajectories  and  different  configurations 
defined  in  the  last  section.  From  these  tables,  it  can  be  seen  that: 

1.  The  P wind  gust  component  does  not  produce  significant  miss 
distances. 

2.  The  lateral  wind  gusts  are  mainly  responsible  for  the  lateral 
errors  for  straight  and  level  flight. 

3.  The  vertical  wind  gusts  are  mainly  responsible  for  the  longi- 
tudinal errors  for  straight  and  level  flight. 

4.  The  lateral  wind  gusts  have  a large  effect  on  the  longitudinal 
error  for  turning  flight  since  a lateral  perturbation  causes  the 
velocity  vector  to  be  perturbed  in  the  vertical  plane  as  well  as 
the  horizontal  plane. 

5.  The  vertical  wind  gusts  have  a significant  effect  on  the  lateral 
error  for  turning  flight  since  a pitching  perturbation  causes  the 
velocity  vector  to  be  perturbed  in  the  horizontal  plane  as  well  as 
the  vertical  plane. 

6.  The  miss  distances  are  significantly  reduced  when  the  aircraft 
Is  flying  in  a descending  turn  because  of  the  reduced  time  of  fall 
of  the  bomb. 
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Table  VII 


Error  Analysis  for  Aircraft  without  Fire  Control  System 

(Pilot  Nulls  1^) 


a.  Oo  - RMS  Value  for  Longitudinal 

Error 

Type  Gust 

Straight  & Level 

Level  Turn 

Descending  Turn 

Lateral 

.3 

18. 

7. 

Vertical 

101. 

54. 

41. 

P 

.2 

7. 

2. 

Composite 

101.0 

59. 

42. 

b.  og  - RMS  Value 

for  Lateral  Error 

y 

Type  Gust 

Straight  & Level 

Level  Turn 

Descending  Turn 

Lateral 

8. 

6. 

3. 

Vertical 

1. 

51. 

15. 

P 

2. 

5. 

1. 

Composite 

9. 

56. 

16. 

c.  Op  - RMS  Value 
®T 

for  Total  Error 

Type  Gust 

Straight  & Level 

Level  Turn 

Descending  Turn 

Lateral 

8. 

19. 

8. 

Vertical 

101. 

64. 

41. 

P 

2. 

9. 

2. 

Composite 

102. 

71. 

42. 

RMS  Values  are  in  Feet. 

C 


r i 

Table  VIII 

Error  Analysis  where  Pilot  Mulls 
a.  Og  - RMS  Value  for  Longitudinal  Error 


Type  Gust 

Straight  & Level 

Level  Turn 

Descending  Turn 

Lateral 

7. 

101. 

40. 

Vertical 

101. 

45. 

39. 

P 

.1 

1. 

.3 

Composite 

101. 

111. 

56. 

b.  0 - RMS  Value 

for  Lateral  1 

Error 

Type  Gust 

Straight  & Level 

Level  Turn 

Descending  Turn 

Lateral 

91. 

57. 

34. 

Vertical 

2. 

39. 

14. 

P 

.3 

.3 

.2 

Composite 

91. 

69. 

37. 

c.  0 - RMS  Value 

for  Total  Error 

^T 

Type  Gust 

Straight  & Level 

Level  Turn 

Descending  Turn 

Lateral 

91. 

114. 

50. 

Vertical 

101. 

53. 

39. 

P 

.4 

1. 

.3 

Composite 

137. 

128. 

64. 
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Table  IX 


Error  Analysis  where  Fire  Control  System 
is  Connected  to  the  Autopilot 

a.  - RMS  Value  for  Longitudinal  Error 


Type  Gust 

Straight  & Level 

Level  Turn 

Descending  Turn 

Lateral 

13. 

127. 

73. 

Vertical 

104. 

45. 

39. 

P 

0. 

4. 

1. 

Composite 

105. 

135. 

83. 

Type  Gust 

b.  a - RMS  Value 

Straight  & Level 

for  Lateral 

Level  Turn 

Error 

Descending  Turn 

Lateral 

75. 

56. 

35. 

Vertical 

1. 

28. 

8. 

P 

0.1 

2. 

1. 

Composite 

75. 

63. 

36. 

Type  Gust 

c.  a - RMS  Value 
®T 

Straight  & Level 

for  Total  Error 

Level  Turn  Descending  Turn 

Lateral 

73. 

137. 

76 

Vertical 

104. 

47. 

39 

P 

0.3 

4. 

1 

Composite 

128. 

147. 

87 

RMS  Values  are  in  Feet. 

C 
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Since  it  is  difficult  to  tell  exactly  which  perturbations  are  responsi- 
ble for  the  miss  distances,  a computer  run  was  made  for  the  level  turn 
case  to  determine  the  significant  perturbations.  Table  X is  a listing 
of  the  rms  values  for  each  perturbed  quantity whi ch  enters  into  the 
error  equation.  For  the  lateral  error,  3 and  i/)  are  the  two  main  contri- 
butors to  the  miss  distance  while  t^  is  the  major  contributor  to  longi- 
tudinal error  with  u contributing  about  20%. 


Table  X 

RMS  Values  for  Error  Parameters  (Level  Turn) 


Variable 

Lateral  Gust 

Vertical  Gust 

P Gust 

°u 

1.67 

2.29 

.01 

Og  (deg) 

.48 

.0 

.0 

(deg) 

.04 

.34 

.0 

% (deg) 

.Z1 

.33 

.0 

Og  (deg) 

.39 

•14 

.01 

0^  (deg) 

2.00 

.14 

.08 

(sec) 

.15 

.05 

.0 

Pilot  Workload 

The  last  area  to  be  considered  to  determine  if  the  pilot  is 
capable  of  manually  accomplishing  the  bomb  run  is  to  determine  the 
Cooper-Harper  rating  for  the  task.  The  Cooper-Harper  rating  is  a 
function  of  performance,  pilot  lead  and  aircraft  response.  Reference 
7 derives  an  equation  which  assigns  a number  to  this  pilot  rating  for 
the  roll  task  being  considered  in  this  study. 


PR  = Ri  (Perf)  + R2  (Tl)  + R3 


(6-10) 


where  PR  = Pilot  Rating 


(Ref  7:16) 


Ri  (Perf) 


R2  (Tl) 


Performance  Contribution 

1.31  0.  Or  "1 
4>c 

Pilot  work  load  due  to  lead  time. 


3.25 


(1-e 


Ratio  of  undamped  natural  frequency  of  numerator 
quadratic  of  ip/6j^  and  undamped  natural  frequency 
of  the  dutch  roll  oscillation. 


= 6.66  1 - -J 

“d 

For  the  F-15  aircraft,  is  equal  to  1.04  for  the  flight  conditions 

used  in  this  report.  For  the  pilot  model  described  in  chapter  5,  the 
pilot  must  generate  a lead  time  of  0.5  seconds.  The  lateral  wind  gust 
defined  in  chapter  4 sets  equal  to  .52  degrees.  The  last  parameter 
required  to  solve  for  the  pilot  rating  is  oa  . These  rms  values  are 
tabulated  in  Table  V where  the  pilot  attempts  to  null  the  command  bank 
angle  41^..  For  the  descending  turn  case,  is  equal  to  .0115  degrees. 
The  pilot  rating  for  this  run  is  determined  to  be  2 where  system  charact- 
eristics are  defined  to  be  good  and  pilot  compensation  is  not  a factor 
to  achieve  the  desired  performance.  For  the  straight  and  level  traject- 
ory where  0.  is  equal  to  5.15°,  the  pilot  rating  is  determined  to  be 

T W 

8.5.  This  rating  indicates  there  are  major  deficiencies  in  the  system 
and  that  considerable  pilot  compensation  is  required  for  aircraft 
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VII.  Conclusions 


If  a pilot  is  to  be  considered  capable  of  releasing  a weapon 
in  a turn,  he  must  be  able  to  act  upon  the  controlled  element  in 
such  a way  so  as  to  substantially  reduce  the  error  signal  without 
causing  an  undue  work  load  upon  himself.  The  transfer  function 
indicates,  that  with  sufficient  gain,  the  simple,  single  axis  pilot 
model  which  was  chosen  in  chapter  5 is  able  to  accomplish  this  roll 
task.  The  pilot  gains  encountered  in  this  study  are  not  excessive 
and  should  not  be  a factor  in  determining  the  pilots  workload.  The 
.5  seconds  lead  time  and  the  .3  seconds  pure  time  delay  are  well  within 
the  normal  response  times  for  a pilot. 

The  pilot  model  is  very  effective  at  nulling  the  errors  in  4)^ 
due  to  wind  gust  inputs.  In  fact,  the  pilot  model  does  a better  job  of 
nulling  the  error  signal  than  the  automated  fire  control  system  which 
ties  in  directly  to  the  flight  controls.  While  reducing  the  error 
signal  the  pilot  tends  to  generate  a larger  perturbed  bank  angle 
than  the  automated  system.  This  indicates  that  the  pilot  inputs  are 
not  as  smooth  or  efficient  as  those  through  the  automated  system. 
However,  it  should  be  noted  that  the  miss  distances,  as  a result  of  the 
larger  perturbations  in  41,  are  small  when  compared  to  the  miss  distances 
caused  by  ip,  b,  and  t^. 

The  miss  distances  for  the  fully  automated  system  and  the  manual 
system  are  comparable.  However,  it  should  be  noted  that  for  the  basic 
system,  where  the  pilot  is  attempting  to  null  4 instead  of  the  miss 
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distances  are  substantially  reduced.  This  indicates  that  the 
wrong  error  signal  is  being  controlled  by  the  fire  control  system. 


Also,  the  miss  distances  for  the  bombs  which  were  released  in  a 
descent  are  much  smaller  than  those  from  a level  flight  condition. 

This  is  due  to  the  fact  that  the  time  of  fall  of  the  bomb  is  shorter 
which  is  directly  proportional  to  the  miss  distance.  It  may  also  be 
assumed  that  even  larger  miss  distances  may  be  expected  for  a bomb 
which  is  released  in  a climb.  One  other  method  that  may  be  used  to 
reduce  the  miss  distance  is  to  reduce  the  release  velocity  which  is 
directly  proportional  to  the  miss  distance. 

The  equations  used  in  this  study  should  be  restricted  to  turning 
flight  where  \ is  not  equal  to  zero.  This  is  the  major  reason  for 
the  large  pilot  rating  for  the  straight  and  level  trajectory  where  the 
system  becomes  too  sensitive.  Otherwise,  the  Cooper-Harper  rating  for 
the  system  indicates  that  it  is  a relatively  easy  task  for  the  pilot. 

If  the  straight  and  level  case  is  disregarded,  the  data  obtained  in 
this  study  indicates  that  the  pilot  is  capable  of  performing  the  task. 
This  finding  is  substantiated  by  the  people  who  have  recently  flown  the 
General  Electric  simulator  which  uses  the  "Firefly"  bombing  law. 


( 
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Appendix  A 


Stability  Derivatives  and  Basic  Aircraft  Data 


The  following  stability  derivatives  are  based  upon  flight 

conditions  of  .6  Mach  at  10,000  feet  altitude  for  the  F-15  aircraft. 

These  values  are  for  body  fixed  axes. 

* = - .8595 

y3 

**  C = .1 


**  C = .32  - .904q,„ 
yr  ° 


Cw  = - .086 


* Cy.  = .1312 

yi5r 


* = - .09741  - .2517ao 


* C,  = - .24 


**  C,  =•  .048  + 6.66ao 


* C»,  - .0373 


* C,  = .04269 

^«DT 

* C,  = .00327  - .00878ac 


* Derivative  obtained  from  Reference  8 for  F-15  aircraft 
**  Derivative  obtained  from  Reference  9 for  YF-15  aircraft 
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* 

C = 
"a 

- .2961 

r = 

a 

- 1.2 

* 

Cm  - 
q 

-3.7 

*4r 

= - .688 

•k 

s- 

.1834  - .467a 

0 

kit 

^np 

' - .03 

k 

Cnr  = 

-.250 

k 

: .00458  - .0416a^ 

k 

= .0229 

k 

^ - .086 

k 

Cl  = 

3.782 

k 

C[)(a) 

= .020  - .0905a  ■ 

Note:  These  derivatives  are  only  valid  for  trim  angles  of  attack  less 
than  8°, 

Values  for  a are  in  radians, 
o 

Stability  derivatives  are  non  dimensional,  therefore,  control 
inputs  must  be  in  radians. 

The  following  data  was  obtained  from  Reference  8 and  represents  a 
typical  configuration  for  the  F-15  aircraft. 

M = 1085,  slugs 

S = 608,  ft^ 
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b = 42.7,  ft 


c = 16.0,  ft 

= 25,274,  slugs-ft^ 

2 

ly  " 155,746,  slugs-ft 
= 175,516.  slugs-ft^ 

Ixz  = - 805,  s1ugs-ft^ 

These  moments  and  products  of  inertia  are  for  body  axes. 
For  the  purposes  of  this  report  the  following  values  were 
considered  to  be  constant: 

P/P^  = .971 

Vt  = 665,  ft/sec  for  .6  Mach 

* O 

Kg  = 2.4  X 10'^  ft  for  low  drag  MK83 
g^  = 32.2,  ft/sec^ 
q = 510.5,  Ib/ft^ 
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Aircraft  Equations  of  Motion 


Nonlinear  Aircraft  Equations 

The  general  force  equations  of  motion  for  the  aircraft  are  based 


m^dV,  /^dV-r  -y  \ ^ 

' = m !•  + ae-v  X Vt  = F 

dt  \ dt  / 


(B-1) 


where  m = Mass  of  the  aircraft 


= Acceleration  in  an  earth  fixed  reference  frame 


Acceleration  in  an  aircraft  fixed  reference  frame 


= Ui  + Vj  + Wk 


-’co  = Angular  velocity  of  the  aircraft  reference  frame  relative 
to  an  earth  fixed  coordinate  system 


= Pi  +Qj  +Rk 


Vj  = Ui  +Vj  +Wk 
f = Fa  t F.  * Ft 


=«  Aerodynamic  forces 


''a,  ^ ^ f^Av  j ^ '"a,  k 


It  should  be  noted  that  i,  j,  and  k are  unit  vectors  for  the  body 
axes  of  the  aircraft.  With  the  aid  of  the  Euler  transformation  matrix 
tJ  (see  Chapter  2),  the  body  forces,  Eg,  may  be  determined  as 
follows: 

Fb-"'[t][o  0 gJT 

= - mg^  sinG  i + mg  cosO  sin4>  j + mg^  cosG  cos4>  k 

where  g^  is  the  acceleration  coefficient  due  to  gravity. 

-► 

Also,  Fj  = Thrust  forces 

= Ft  i + Ft  j + Ft  k 
•x  Ty  'J  Tz 

X V = (QW  - RV)  i + (RU  - PW)  j + (PV  - QU)  k 


Using  these  definitions,  it  is  possible  to  write  the  general  force 
equations  of  motion  as: 

f'Av  •'t 

U + QW  - RV  = - g_  SinG  + _ + 

mm  (B-2) 

. F^  Fj 

V + RU  - PW  = g_  COSG  sin$  + — ^ + — L 

^ mm  (B-3) 

. F/\  Fj 

W + PV  - QU  = g_  COSG  COSO  + 1 + L 

""  m m (B-4) 

(Ref  10:2.24) 


The  rate  of  change  of  the  angular  momentum,  as  seen  by  an 
observer  in  inertial  space,  is  equal  to  vector  sum  of  the  external 
moments.  Mathematically  this  can  be  written  as 


edH 

dt 


— + ae;*;  X ff  = M 

J..  w 


(B-5) 


) 
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where  H = Angular  nromentum  vector 


XX 

~^xy 

-hz 

yz 

^yy 

•^yz 

zx 

~^zy 

^z 

P 

Q 

R 


= Rate  of  change  of  angular  momentum  in  the  aircraft 

dt 

reference  frame 

= - 'xz«)  ' ♦ ('yy«  i * <'zz«  - 'zx^l  k 

where  the  xz  plane  is  a plane  of  symmetry  and  = I = “ 0- 

yx  y^  ^y 

X it  . [-  PQ  4 Qr]  f 4 PR  4 

l„  (p"  - r")]  i 4[l„  QR4  (lyj,  - pqj  £ 

= M;^  + ^ 

Ma  = l-A  ^ ^ ”a  J ^ Na  k 

i + j > Nj  ic 

Combining  the  above  equations  yields: 

‘xx  ^ ■ 'xz  ' - 'xz  '■'>  " f'zz  - 'yy'  W ’ '■A  * '■T  ("-«) 

'yy  « * f'xx  - 'zz)  >’«  * 'xz  f'  - * "a  * H,  <8-') 

'zz  “ - 'xz  P * <'yy  - ‘xx)  '■')  ♦ 'xz  w ■ "a  * "T  I'-®) 

(Ref  10:2.10) 
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These  six  equations  of  nwtion  (equations  B-2  thru  B-4  and  B-6 
thru  B-8)  contain  eight  unknowns  which  are:  U,  V,  W,  P,  Q,  R,  0,  <P. 


To  correct  this  problem,  the  three  kinematic  equations  may  be 
derived  with  the  aid  of  the  Euler  transformation  matrices  (see 
Chapter  2). 


1 0 -sin© 

* 

0 cos*  cos0sin* 

9 

0 -sin*  COS0COS* 

y 

Solving  the  above  linear  set  of  equations  for  <t,  0,  and  i'  yields: 


$ = P + Q sin4>  tans  + R cos*  tane  (B-9) 

0 = Q cos4>  - R sin*  (B-10) 

y = (Q  sin*  + R cos*)/cos0  (B-11) 

(Ref  10:2.22) 

In  the  development  of  the  previous  equations,  the  following 
assumptions  apply: 

The  earth  fixed  coordinates  are  considered  to  be  inertial 
coordinates. 

The  aircraft  coordinate  system  is  body  fixed  with  the  origin  at 
the  center  of  gravity. 

The  mass  of  the  aircraft  is  constant. 

The  angular  momentum  due  to  spinning  rotors  has  been  neglected. 

The  xz  plane  is  a plane  of  symmetry. 

The  aircraft  is  a rigid  body.  ' 
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In  order  to  solve  these  nine  aircraft  equations  of  nation,  it  is 

necessary  that  U,V,W,P,Q,R,o,  and  <ii  be  known.  From  the 
O 0 O 0 0 0 o o 

geometry  equations  (see  Chapter  two),  values  for  a^,  and  were 

determined  by  solving  for  the  required  trajectory  of  the  aircraft. 

Knowing  a^,  and  can  be  obtained  from 

U = Vt  cos  a (B-12) 

o lo  o 

W = V-r  sin  a (B-13) 

0 Tq  o 

Also  for  coordinated  flight  with  zero  sideslip,  V = F.  =0. 

y 

Equations  B-2  through  B-4  and  B-6  through  B-10  may  be  used  to  solve  the 

remaining  eight  unknov/ns  (P  » Q^.  N^)  for  steady 

X z 

state  flight  conditions.  P , Q , R can  be  obtained  from  equations 
B-3,  B-9  and  B-10  as  follows: 


Q 


o 


P 

0 

R 

0 


= R tan  (t 
0 0 


= - R^  (tan(}i^sin(ii^tano^  + cos(j>^tan0^) 

g sin*  cose 
C *^0  o 


U + W (tan*  sin*  tane  + cos*  tane  ) 
0 0 ^0  ^0  o o o 


(B-14) 

(B-15) 

(B-15) 


Linear  Aircraft  Equations 

These  nine,  coupled,  nonlinear,  differential  equations  describe 
the  general  motion  of  the  aircraft.  Since  this  study  is  interested 
in  the  aircraft  dynamic  frequency  response  about  some  equilibrium 
state,  the  perturbation  technique  is  appropriate.  Using  this  approach, 
the  equations  may  be  linearized  and  Laplace  transformations  used  to 
study  the  frequency  response.  This  can  be  accomplished  by: 

(1)  Considering  the  eight  state  variables  (U,  V,  W,  P,  Q,  R, 
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0,  ■>)  as  being  composed  of  some  equilibrium  (steady  state) 
value  and  some  perturbation  value.  Those  variables  can  then 
be  written  as: 


u = 

u + 

0 

u 

P = 

P + 
0 

P 

0 = 

Q 

0 

+ e 

V = 

V + 

0 

V 

Q = 

q 

4>  = 

<P 

0 

+ 4, 

w = 

w + 

w 

R = 

R + 

r 

The  thrust  and  aerodynamic  forces  and  moments  may  be  expressed  as; 

Ft  = Ft  + f7 


^Ax  ^ ^Ax  '^Tx  "^Tx  'Tx 

o o 


S ° So  S ''Ty  ■ s„  * s 

Z 

^T  " ^To  *T 


f^Az  = ^Az„  ^A, 


•-A  = >-«  ^A 


^A  = \ + "’A 


h = \ ^ "A 


Mt  = 


Nj  = Nt  + nj 
* /^ 


Where  the  steady  state  values  are  denoted  by  the  subscript  (o) 
and  the  perturbed  state  quantity  is  represented  by  lower  case 
symbol . 

(2)  Defining  the  perturbation  angles  <()  and  e such  that: 

cose  = 1 coS(j>  = 1 

sine  « e sin(>  = 

In  general,  this  relationship  holds  for  angles  less  than  15°. 
(Ref  10:2.33) 
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'.-iv-.'" 

i-  t 


(3)  Assuming  that  all  products  and  cross  products  of  the  pertur- 
bation variables  may  be  neglected  when  compared  to  the  magni- 
tude of  the  perturbations  themselves  or  by  using  the  linear  terms 
of  a Taylor  series  expansion  about  the  nominal  conditions. 

(4)  Noting  that  the  steady  state  equations  of  motion  and  kine- 
matic equations  which  are  embedded  in  the  perturbed  state  equa- 
tions have  already  been  satisfied  and  therefore  can  be  eliminated 


(Ref  10:2.33). 

With  these  basic  assumptions,  the  linearized  force  equations  are: 


u-Vr-Rv  + Wq  + Qw  = -g  ocose  + — - + 

O 0 0 0 ^c  o_  _ 

m m 


^ (B-17) 


v + Ur  + Ru-Wp-Pw  =-g^  esine  sin^  + 
0 0 O 0 C 0 0 


% . ^ 


g ^cose  cos(j)  + + 


m m 


(B-18) 


w-Uq-Qu  + Vp  + Pv  = -g^  esine  cosa  - 
0 o o'^  o ^C  o o 


g^.  (j.CDS0^sin(>^ 


(B-19) 
(Ref  10:2.33) 


For  this  analysis,  the  perturbations  due  to  thrust  are  assumed  to  be 
small  and  have  been  neglected.  Also,  for  normal  flight  conditions,  the 
aircraft  is  assumed  to  be  flying  a trajectory  that  will  result  in  a 
reliable  release  without  any  inputs  to  the  controls.  This  assumes 
the  aircraft  is  in  coordinated  flight  where  V is  zero. 

^Ay 

The  value  -jjp  may  be  determined  from  the  following  equation: 

F.  » C 5 S (B-20) 
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where,  = C|_  (a)  sina  - (u)  cosa  - 3.702a  sina  - (.020-. 0905a  + 

(Ref  8 : 19,20) 


2.325a^)  COSa 


This  equation  may  be  linearized  by  using  the  Taylor's  series  expansion 
limited  to  linear  terms  where: 


\^\  + 


H 3a 


1 0 


3Fa 

where, — - 
a/u 

lU 


q S 3C„  C„  S 3q 
1.+  ^ 


3(  u 

lU 


lu, 


2 q S C, 


3*^ Ax  ^ ^ 3^X1 


3a 


3a 


o o 

f A 2 q S u 

x = (3.782a  sina  - .020  COSa  + .0905a.  COSa.- 

m m U 

O 

2 ^ 

2.325a  cosa  ) + (3.782  sina  + 3.782a^  cosa^  + 

rt  » • r\  O O 


m 


.020  sina  + .0905  cosa„  - .0905a  sina  - 4.65a  coSa  + 

O O V O o o 


2.325a  sina  ) 


(B-21) 


Combining  like  terms  and  using  the  flight  conditions  referenced  in 
fA 

appendix  A,  — ~ may  be  written  as: 
m 


— = u r.86  (3.782a  sina  - .020  cosa  + . 
L 0 0 0 

m 


0905a  COSa  - 
o o 
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2.325a  ^ cosa  )1  + a r286  (3.802  sina  - .8680a  cosa  ♦ 

0 o J L 0 0 0 

.0905  cosa  - .0905a  sina  + 2.325a  ^ sina  )!  (B-22) 

O 0 0 o o J ' 


The  value  for — =-  may  be  determined  using  the  same  method, 
m 

q S 


(B-23) 


-2. 


where,  ~ \ o - Cq  (i)  sin  a 

* - 3.782i  cosa  - (.020  - .0905a  + 2.325a")  sinS 
Using  the  Taylor's  series  expansion. 


% . 

" IS) 
0 

o 

u q S aC^ 

= 2 q S C,  — + — 

^ U aa 

0 


= u r.86(-3.782a  cosa„  - .020  sina  + .0905a  sina„  - 

m L 0 0 0 0 0 

2.325a^^  o[286(-3.802  cosa^  - .BGOa^^sinag  + 

.0905a  coSa  .0905  sina.  ~ 2.325a  coSa  )1  (B-24) 

0 0 o 0 0 J 


y may  be  determined  from: 

m 


m m 


C 


C pb  C , rb 

«#  • \t 

2U 


‘y  = __  / >6  6 + yp  — + 'yp  + c 6^^  + 

2U  ^®DT 
o 


(B-25) 
(Ref  10:4.113) 


286 


. [-  .85958  + .00321p  + (.0103  - .029ao)r- 


.086  6qj  + .1312  6pj 


(B-26) 
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Combining  equations  B-17  and  B-22,  B-18  and  B-24,  and  B-19  and  B-26 
and  also  collecting  like  terms  yields: 

uT.Se  (3,782a  sina  - .020  COSa  + .0905a  COSa  - 2.325a  ^ COSa  ) - 
L 00  o oo  00 

sj  + 6(665  R^)  + 0^286(3.802  sina^  ■ .8680^  cosa^  + ,0905  cosa^  • 

.0905a  sina  + 2.325a  ^ sina  ) - 665  Q 1 - q (W  ) - 
0 o o 0 ' qj  ^ o ' 

0 (32.2  cosej  = 0 (B-27) 

- u (R^)  - 6(245.8  + 665s)  + a(665P^)  + p (W^  + .918)  + r (2.946  - 
8.249ao  " 0^)  - 9 (32.2  sino^  sin({i^)  + <ti(32.2  cose^^  coscji^)  - 

(24.6)  + 6^  (37.52)  - 0 (B-28) 

u To  + .86  (-  3.782a  cosa  - .020  sina  + .0905a  si  no  - 
L o 00  0 00 

2.325a  ^ sina  )!  - 6(665P  ) + ar286  (-  3.802  cosa  - .868a  sina  + 

0 o'J'o'L'  0 00 

.0905a  cosa  + .0905  sina  - 2.325a  ^ coSa  ) - 665sl  + q (U  ) - 

00  0 0 o J o 

0 (32.2  sine  C0S(t>  ) - (j.  (32.2  cose  sin*  ) = 0 (8-29) 

0 O O 0 

It  should  be  noted  that  in  the  proceeding  equations  is  equal  to 
zero  and  s represents  the  Laplace  transform  variable. 

The  linearized  moment  equations  may  be  written  as: 

^ P - ^xz  ^ - ^xz  (^0  ^ ^ ^ ■ ^y^ 

(8-30) 

ly  q + dx  - ^z)  ^ ^xz  (2  P^P  - 2 R^r)  = m^  + mj  (8-31) 
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I,  r - Ix,  p ^ (ly  - I,)  (P,  c,  ^ Q,  P)  ^ Ixz  (Q,  r q) 


= (B-32) 

(Ref  10  ;2.33) 

The  same  assumptions  apply  here  as  in  the  force  equations,  that 
is,  the  perturbations  due  to  thrust  are  assumed  to  be  small  and  = 0. 
The  perturbed  rolling  moment  may  be  written  as: 


= q S b(c,^  6 + C 


pb  rb 


~ + C„  Co . 6a  Co  6r.T 

*-6a  a Jl5Qy  ^DT 


ilp  2U  2U 

o o 


'Ur  ^r) 


(B-33) 


(Ref  10  :4.113) 

Again,  using  the  flight  conditions  from  appendix  A: 

= 13.25  X 10^  [(-  .09741  - .2517a^)  3 - .0077  p + (.00154  + 

.ZlSSa^)  r + .0373  + .0427  6^.^+  (.00327  - .C0878a^)  6^  (8-34) 

The  perturbed  pitching  moment  may  be  determined  as  follows: 


_ _ u aC  qc 

m A = q S c (Cm  “ + C a + C^-  + C„  + 

rt  "’u  u "la  \ 2U  <1  2U 

O 0 0 


(B-35) 
(Pxef  10:4.113) 


where  C-,  is  negligible  for  mach  numbers  below  the  transonic  range. 


m^  = 4.96  X 10  (-  .2961(1  - .0144(i  - .0445q  - .688  6^j) 


(B-36) 


The  perturbed  yawing  moment  is: 

pb  rb 

• C'. 

2U 


n*  = q S b 


( ^nr, 

\ 3 P 2U 


^ ^nsoT  ^UT  ^ 

0 a 


(B-37) 
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= 13.25  X 10*^  [(.1834  - .467«^)e  - .000963  p - .0080  r + 

(.00458  - .0416a^)  6^  + .0229  6^.^  _ ggg  6pJ  (B-38) 

After  combining  the  terms  for  each  variable,  and  dividing  through  by 
6 

10  , the  three  angular  nomentuin  equations  are: 

-6  (1.29  + 3.335a  ) - p (.102  + .0253s  + .0008  Q ) - q (.0008P  + 
o 0 0 

.0197  R ) + r (.0204  + 2.833a  - .0008s  - .0197  Q ) + 6»  (.494)  + 
o 0 o A 

(.5656)  + 6,^  (.0433  - .116a^)  = 0 (B-39) 


-a  (1.47  + .0717s)  + p (.0016  + .150  R^)  - q (.221  +. 1557s)  + 


r (.150  P„  - .0016  R ) - (3.416)  = 0 

o 0 ti  I 


(B-40) 


6 (2.43  - 6.187a^)  - p (.01276  + .0008s  + .130  Q^)  + q (-.130  P^  + 
.0008  R ) + r (-  .1063  - .1755s  + .0008  Q ) + (.0607  - 

0 0 M 


551a J + 6qj  (.303)  - 6p  (1.14)  = 0 


The  three  linearized  kinematic  equations  are: 


sin(ji  coS(j 

q r - ij/  s +9 

cose  cose 


tane 

?(Q  sin*  + R cos*  ) 

COSO  o 0 0 o 


(B-41) 


/ COS* 

^ ( Qo ^ - R, 

V COS0„ 


sin(j> 


0 


cose 


= 0 


q (cos<|)  ) - r (sin(|>)  -os  - iji  (Q  sin^  + R coS(}>  ) = 0 


O 0 


(B-42) 


(B-43) 
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t 

I 


1 


p + q (sin<{i  tano  ) + r (cosi}i 

0 0 o 


cos^o 


2^  I + Ij)  (Q^  cos4>g  tano^  - tano^  - s)  = 0 (B-44) 

° cos  Oq/  ° 00 

The  nine  linearized  aircraft  equations  are  summarized  in  Table  XI. 


Table  XI 


ao>ao«6  aaa 


UNCLASSIFIED 


AIR  FORCE  INST  OF  TECH  MRI4HT-PATTERS0N  AFB  OHIO  SCH-- ETC  F/9  15/ 7 
A FEASIBILITY  STUDY  OF  A MANUAL  BOMB  RbLEASE  WHILE  IN  A TURN.(U) 

DEC  77  U 0 WALTON 

AFIT/BAE/AA/77D-17  NL 


2 -78 


Table  XI  (Continued) 


cos  6/  \ cos 


Appendix  C 


Flight  Cont roi  Equations 

Since  this  study  is  based  upon  the  frequency  response  of  the 
aircraft,  linearized  flight  control  equations  must  be  included  in  the 
system  model.  The  F-15  flight  control  system  is  basically  linear  with 
only  a few  appropriately  placed  limiters,  deadbands  and  gain  schedulers. 
Because  this  study  is  an  examination  of  the  dynamic  response  about  some 
nominal  flight  condition,  most  of  these  nonlinearities  can  be  either 
avoided,  neglected  or  linearized. 

The  actuators  for  the  control  surfaces  are  depicted  in  figures 
11  through  14  (Ref  11  :25,26).  Since  these  actuators  are  capable  of 
responses  of  up  to  20  radians  per  second,  they  have  not  been  included 
in  this  analysis.  The  reason  for  this  is  the  primary  frequency  range 
of  interest  is  below  ten  radians  per  second.  The  simplified  block 
diagrams  for  the  pitch,  roll  and  yaw  axis  are  presented  in  figures  15  . 
16  , and  17  (Ref  11:29-32).  It  should  be  noted  that  all  the  inputs 

are  in  pounds  force  and  all  the  control  displacements  are  in  radians. 

F-15  Flight  Controls 

With  these  assumptions  and  clarifications,  the  primary  pitch 
axis  control  laws  may  be  written  as: 
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Input 


HT 


Fig.  13.  Horizontal  Tail  Actuator 


Roll  Mech 
Input 


Roll  Elect 
Input 


Fig.  14.  Differential  Tail  Actuator 


Pitch  Elec  Input  =< 


f 

( 1 \ 

( ^ \ 

p 

- 

'"es  ' 

\2.73l 

u 

+ 

i + q (12.42)  ( — ] + 
mg^  \3s+l/ 


s q (.68) 


/30+3s\] 
\30+s  /J 


Pitch  Mech  Input  = 


-F„(.l)(2)  + q s (.5) 


s/  \57.3, 


3.75 


fe)  (i) 


(12.42)  ( — )* 
mg  \3s+l/ 


L ■'c 


s q (.68) 


30+3s 


30+s 


1.5 


(C-1) 


} 1.5[l+-l(.2)( F^3(.1)(.6464)(2) 


1 

57.3 

(C-2) 


From  Fig  13.  it  can  be  seen  that: 
Pitch  Elec  Input  + Pitch  Mech  Input  = 


(C-3) 


^HT  " ■''eS 


. 0023s .048s^  + .0526s  + .0105 
s^  + 5s^ 


F^^  /. 0786s ^ +. 8816s +.835s  +.2356 


mg 


3 2 

s + 30s 


,1603s‘*  +4. 798s  ^ +36il4s^  +40.24s+8.94 
3s^  + 91s^  + 30s 


(C-4) 


where = u F. 0267 (-3. 782a  cosa  - .020sina^  + .0905o  sina  - 

I ' 0 0 O 0 0 

2.325o  ^ sina  al8.885(-3.802cosa  - .868«  sina  + 

0 o'J  L 0 0 0 

.GOOSagCOSOg  + .0905sinaQ  - 2.325«^^  (C-5) 

From  Fig.  16,  the  roll  axis  flight  control  equations  are: 
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Fas  "" 


(C-6) 


1 \ / 1 


AS 


1+.3S/  \57.3) 


(.0625)  = 6 


DT 


ELEC 


(C-7) 


Fas  (-24)  (10) 


1 \/.3 


= 6 


.57.3/ \ 2 


DT 


MECH 


(C-8) 


The  two  differential  tail  inputs  (Snj  and  6r,j  ) may  be  summed 

ELEC  MECH 

together  to  obtain  the  differential  tail  displacement  (see  Fig.  14  ). 


'^DT  ^ “^DT  “ ^ '^DT 

‘^'elec  ^‘mech 


(C-9) 


AS 


^0019s  + .0074 
L 3s  + 1 


- .0625p  = 2 5 


DT 


(C-10) 


The  primary  yaw  axis  flight  control  equation  can  be  obtained  from 
Fig.  17  . 


^'^\57.3AiO+s 


—y(. 00537)  + rf— ] - (aP  +pct  )f— 


m \2s+l/  ° ° \s+l; 


S+.5 


1 


(.84)  + Fpp  (.5)  = (C-11) 

s / 57.3 


•^A  r 

where  1=  286  -.8595e  + .00321p  + (.0103  - .029a  ) r - 

m L 0 


.086  Sqj  + .1312  6pJ  (C-12) 


.0087s  + .174\  /1.109s  + .555 


Rpl 


s + 10 


= B 


r.84P  s + .42 PA 

o o \ 

+ a| + 

S + 1 / 
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’(.84a  - .0041)  + (.42a  - .0062)s  - .002l\ 

' 2 2 + 


+ s 


(.15a  -3.4)s^  +(.22a  -1.76)s^  ^■(•IK  -.04)s+.019a„  -.0066' 
r I 2 0 ^ ^ 0 o , ^ 

4s  + 4s  + s 


.1107S+.0553 


.8307s -.0845 


^DT 


+ 6r 


(C-13) 


Firefly  Flight  Controls 

The  "Firefly"  automated  system  has  been  studied,  in  order  to  have 
a baseline  with  which  to  compare  the  manually  operated  system.  Figures 
18  , 19,  and  20  depict  the  F-15  flight  control  lav;s  whicli  have  been 

modified  for  the  Firefly  configuration.  Referring  to  Fig.  18  , the 
augmented  pitch  axis  flight  control  equation  may  be  written  as: 


Pitch  Elect  Input  = 


\5+s/  \30+s  / 


(C-14) 


Pitch  Mech  Input  = | 


-Fp  (.1)(2)  -—  + q s (.5) 

'"9c 


3.75 


ES 


3.73/V  5+s 


^ + q (12.42)  f— ) + 
l'"9c 


(C-2) 


s q (.68) 


’30+3s 

> 

i 

/ l\  1.5 

1.5  1+-  (.2)  Fp.(.l)(.6464)(2) 

J 

_30+s  _ 

> 

> \ s/  s 

57.3 
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By  summing  the  electrical  input  and  the  mechanical  input,  the  pitch 
axis  equation  for  the  flight  control  is: 


^HT  " ‘^ES 


,0023s^  + .0136s^  + .0196s  + .0105' 
__  - 

5s  + s 


Fa  /.0166s  + .049s  + .2356’ 


mg. 


2 3 

30s  + s 


.158s^  + 3.06s^  + 13.82s  + 9.94’ 

*4  I 3 2 

3s  + 91s  + 30s 


(C-IS) 


For  the  aileron  input  there  is  no  change  from  the  basic  F-15  flight 
control  system.  Therefore: 


F,s  = 23.87 


{C-16) 


However,  the  output  from  the  fire  control  equation,  p^.,  is  fed 
into  the  roll  axis  of  the  differential  tail.  Figure  19  shows  the  p^, 
input  plus  the  slight  change  in  flight  control  configuration.  This 
diagram  leads  to  the  following  equations: 


Roll  Elect  Input  = 


Roll  Mech  Input  = Fas(.24)(10) 


(.365) 


(C-17) 


(C-18) 


Recalling  that: 

Roll  Mech  Input  + Roll  Elec  Input  = 2 <5^^  (C-19) 


88 


.0019s  ' +.0127s4-.0128\  /.365s>.73  1.025s*3.65 


AS 


s + .3s 


V 


5s  + s" 


=2  5, 


(C-20) 


The  augmented  yaw  axis  flight  control  equation  can  bo  obtained  from 
Fig.  20. 


.0087s+.17^  L.035S+.5182' 


RP 


s + 10 


+ aP 


,78s+.39' 


s + li 


(.78a  -.0039)s  +(.39a  -.0058)s-.0019' 
0 0 

s^  + s 


/(.14a  -3.41)s^  +(.21a  -1. 755)s^  +(. 105a  -.0376)s+(.0175a„  -.0063)^ 

I 0 0 0 o ' 


4s ^ + 4s^  + s 


+ 6 


DTl 


103s+.0517\  /.776s-. 0789 

" ^R 


(C-21) 


( 
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Fig.  19.  Augmented  Roll  Axis  Control  Laws  (Firefly) 


Fig.  20.  Augmented  Yaw  Axis  Control  Laws  (Firefly) 


Appendix  D 


Computer  Model 

The  computer  program  "FRQRSP"  is  a general  purpose  frequenc;y 
response  program  which  can  be  used  to  evaluate  a system  model  and  plot 
amplitude  and  phase  response  of  a linear  transfer  function  over  a 
specified  frequency  range.  It  is  also  capable  of  solving  a set  of 
simultaneous  equations  by  the  Gauss  Jordan  elimination  method 
(Ref  12:17).  Additionally,  the  program  is  set  up  to  perform  spectrum 
analysis.  The  evaluation  of  the  system  is  accomplished  by  first 
defining  a function  G(s)  where  s is  the  Laplace  transform  variable. 

The  transfer  functions  in  this  report  have  been  defined  as: 

1.  The  open  loop  response  of  the  system  {<^  and  to  a lateral 
stick  input. 

2.  The  open  loop  response  of  the  system  (<j)  and  to  wind  gust 
inputs. 

3.  The  closed  loop  response  of  the  system  (<j(  and  to  wind 
gust  inputs  (pilot  in  the  loop). 

4.  The  closed  loop  response  of  the  system  to  wind  gust 
inputs  (automatic  fire  control  system). 

In  order  to  model  the  system  for  the  computer  program,  it  is  necess 
ary  to  set  up  the  following  complex  matrix  equation: 

Wl"]  •[“] 
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where  [Aj  - Represents  the  system  components  which  consist  of  the 
aircraft,  the  flight  controls,  and  the  fire  control  law 
(Figure  21  provides  a listing  of  the  matrix  elements  while 
figure  22  depicts  the  position  of  the  non  zero  elements  in 
the  matrix). 

- The  unknown  vector  representing  the  system  variables. 

- [u  6 a p q r 0 «A  «DT  «R  V ^ ‘'p„  Pc]^ 

^ B j - The  vector  which  represents  the  external  inputs  to  the 
system  such  as  wind  gusts  and  lateral  stick  inputs  (see 
figures  23  and  24  ). 

The  abbreviations  used  to  define  the  A and  B matrices  in  the  computer 
program  (Fig.  21  , 23  and  24  ) are  listed  below  in  terms  of  the  sym- 
bols used  in  this  report. 

Fixed  inputs: 

ALAM0  = X 

o 

DK  = K|3 
E = b 

GC=  g, 

H = h 

K0  » K 

Ranges : 

G0  = G, 

RP0  = Rp 

r 0 

RR0  = R^^ 

RTV  = Rt 

' V 

R2  = R2 


PI  = 

TT 

RHO  = 

p/p 

RT0  = 

"To 

TF0  = 

VT0  = 
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Velocities: 


U0  = u 

0 

P0  = 

P 

0 

V0  = V 

0 

Q0  = 

Qo 

W0  = w 

0 

R0  = 

R 

0 

Bombing  Law  Components: 

CPW0  = Cp 

RDZ0  = 

DTN0  = D 

' '’o 

WOMEGA0  = 

% 

PC0  = P 

Co 

0 

Angles: 

A0  = a, 

F0  = 

B0  = 3^ 

S0  = 

K 

G0  = 

T0  = 

0 

0 

Trigonometric  functions: 

SA0  = sinog 

CA0  = 

cosa 

0 

SB0  = sin3^ 

CB0  = 

coss 

0 

SF0  = sin* 

0 

CF0  = 

COSA 

0 

SG0  = sinv 

0 

CG0  = 

COSv 

0 

SL0  = sinx 

0 

CL0  = 

cosx 

0 

SS0  =«  sin>p 

0 

CS0  = 

COSlJ;^ 

ST0  = sine 

0 

CT0  = 

cose 

0 

To  determine  the  system  variables  in 

r],  the  A matrix 

must  be  inver 

ted.  This  procedure  uses  pivoting  in  order  to  maintain  an  accurate 


solution.  Note  that  the  ICOL  statements  following  each  element  place 
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1 


the  element  in  the  correct  colunui  for  that  row. 

The  parameter  P(i)  is  used  to  correlate  the  inputs  with  the 
correct  data.  For  this  program,  the  P parameters  are  defined  as: 

P(l)  - Pilot  gain. 

P(2)  - Longitudinal  input  (not  used). 

P(3)  - Lateral  input. 

P(4)  - Horizontal  wind  gust  input. 

P(5)  - Vertical  wind  gust  input. 

P(6)  - P wind  gust  input. 

P(7)  - System  gain. 

To  obtain  the  frequency  response,  the  user  supplied  complex 
function  G(s)  is  evaluated  for  a region  of  frequencies  defined  by 
the  input  data.  The  frequency  is  reported  in  both  radians  per  second 
and  hertz  in  tabular  form,  however,  the  Bode  plots  are  in  hertz.  The 
amplitude  of  G(s)  is  reported  as  magnitude  and  decibels  where  decibels 
are  defined  to  be  20  logig|G(s)|.  The  phase  of  G(s)  is  reported  in 
degrees  as  the  arctangent  of  the  imaginary  part  divided  by  the  real 
part  (Ref  12  :4). 

For  the  PSD  analysis,  the  user  must  again  supply  the  complex 
transfer  function  G(s)  where  the  input  PSD  is  a real  function  in  terms 
of  radians  per  second  or  u.  The  results  of  the  PSD  analysis  consists 
of  a tabulation  of  the  input  and  output  PSD  in  both  magnitude  and 
decibels  where  the  dB  values  are  defined  as  10  logio|PSD|.  The  program 
also  tabulates  and  plots  the  cummulative  input  and  output  power  which 
is  the  square  root  of  the  input  and  output  PSD  magnitude  values  inte- 
grated from  the  starting  frequency  to  the  current  frequency.  When  the 
program  completes  the  evaluation  of  the  input  and  output  PSD  values 
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over  the  specified  range,  the  final  cutmnulative  power  values  are 
reported  as  the  root  mean  square  (RMS)  values  (Ref  12  :5). 


BEST  AVAILA8IE  COPY 


n--IMF  1 IN  riJLI.  MATRIX  FORM 

n:  A(?1,  ?Q) 

ro  i"LFx  A,<;  . 

$ rs  THT  LAPLAOe  TRAMSFCKM  VARIABLE 

A ( I , 1)  =0  .B5»  ( 5 .ri’^AO’S;  C-.  :20»CA0*  .0'rJ5»A'J*CAQ-2. 32?*  Au'‘*?»CA0)  -S 

r.f ':oL<  1,  i)  =1 


TCOl.  II  ,2)  =2 


3)  =?‘3?.  * (.5.  A32*SAC  - . '’0-J»Arj»CAa*  . 09a5*CA0-.0B3‘;*4Q*^ A 0*?.  725* 
2*'3A0)-(Sf)?.*no  5 irnL(l,3)=3 

‘♦)=-WC  .3  rCOL(l,‘)=Si 

5)  =-32.2>-CT0  5 I'.OL  ( 1 , 5)  =d 

1)  =-RC  » tCOL  (?,  n = i 

2)  (2't''-.  'stCCL<2,2)  =2 

3) =3?5.*'’3  STOOL  (2,3)  =3 

'O  =W2f  ICOL  (2,i»)  ='> 

?)  =?.9i.=.-^.2T.*A0-UC  i IC0L(2,5)=6 

=-72.2»>r  n^OFO  S rCOL(2,6)=3 

=32.2' CM ‘C^O  S ICaL(2,2)=q 

3)  =-2'+.6  JICOL  (2,3)  =12 

0)  =37.52  i TOOL (2,0) =13 

1)  =r)0».?F  • (-3.762*  Ar*CAO-.fl2a*SAO  f.C9:0*AC»SA2-2.  325*A0  **2*SA0. 

« ISO'.  ( .3,  1)  =1 

2)  =-?CP. '•■■’0  $TCCL(3,2)=2 

3)  =2  3f>.*  (-3,3J2*CAC-.,30.“‘AO*SAO*.OOC?*AO*CAO*-.  j305*SAa-2.325* 

2*CA0)-f>6;.'S  S ICOL(3,3)=3 

•^)='J0  S iroi.  (3, ;.)  =5 

5) =-22.2‘Oro*CFO  5 ICCL(3,0)=3 

•S)  =-32.2*Cr  j‘3FC  3 IC0L(3,F)=0 

1)  =-l . 21-’.  •'3F*A  3 J ICCL(4,1)=2 

2)  =-.lC2-.32;3*S-.C.:oa*OJ  I TOOL  (!,,2)='f 

3)  =-.0  00  3 0197  SICOL  ('>,3)  =0 

<♦)  =.  0?;-*2. 3 33*  AJ- .1  C 0 3*3-.',  107*  no  3 irOL(^,»)=5 
5)  =.'»0',  l'ICOL('t,‘')  =11 

',)  =.5‘i'6  f IC0L(',,&)=12 

7)  =.G',33-.  lir,*AQ  5 ICCL  ( '* , 7 ) =1 3 


1) =-l.i»'-.1717*S 

2)  =.  3C1‘, *■’)*. I'jO’RO 

3)  =-.22l-.lF0,-*S 

',)  =.15C'P0-.0.)17,*FC 
5)  =-3.VlC, 


3 ICOL(0,l)=3 
5ICOL(5,2)='» 

$ ICOL(!i,3)=0 
31COL  (5,1,)  =6 

3 ICOL(5,O)=10 


1)  =2.1,3-?. 167*A0  f IC0L(6,l)  = 2 

2)  =-.C127',-.l  laS’E-.l ’i)*nj  $ ICOLI6,2)='» 

5)  =-.i3:*>o,.:ooi*Rr  yicol(o,3) =5 
^)=-.1007-.120?*S*.r:  jjvio  J tC0L(6,i.)=6 

5)  =.0?0/-.  OOl’Afl  3 IC0L(?,5)=H 

0)  =.703  J rC0Llo,0)=12 

7)  =-l.l«,  3 lCCL(r>,7>=13 


Fig.  21.  Elements  of  the  A Matrix 
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BEST  AVAIUBLE  COPY 


4 f 1>  ^S'-O/CTJ 
ftt/.c’)  -Ct-'Q/CTO 
M ' , ,T)  =-S 


I rcoi  ( 1) 

J ICrL(/,2)=6 

JiroL  (/,.n=7 

ftf’.'t)  =(rA'uro) /c:o)*(0(:*‘ipo+-{o*:fo)  s ^co^.(7,'^)=^ 
/.(’.j)  =10''CrC/cro-P.O*SF!./CTC  3 rCOL(7,5)=9 


Af<,l)=:PO  7 TCOKJ.Ds*} 

a(),;?)=-SFO  7 TCnL(?,7)=6 

4{i,,7)=-3  3 rCOL(i»,3)=a 

4(i,;»)  =-OC‘SF:ff’u‘CFO  $ K0L(6,t.)=0 

ac*,’)  =1.  » IC0L(9,1)='> 

a(  »,’)  =3'^(1'TA-|(T0)  SrCOL  (0,2)=!5 

An,  5)  =:F0»’'A-HTa)  J rC0L(9,3)=fi 

A f 1 , .)  =3  1 ►SFC/niO  * *?V'’C’ CFu/CTO*  *2  5 TCOL(9,'»)  = 3 

‘\(  1,  •'.)  rlJ'CFC-' TA'KTO) -PC  ♦SFQ'^TAm  TO) -3  S ir0L(0,5)-0 


A(n,l)=  .02f)’*(-3.732»r.C»Cf.  3-.02O*3A0«-.riOJ5-»A0-'3AJ-2.32‘3*Aij»»2» 

C -•  too.  (10,1)  = 1 

A(lC  ,2)  ■:  '’.3')'‘(-3.e]2*rAC-.a‘5)*'AC*OA.jnC135»Aa»CAO*-.0  9a5*OAO- 
0 7. 52i»  Aa*»2»  r\  J ) ‘ ( .C7S6'S»*3!-.  1016*?**2f  .).'3‘3*S).230i»  / ( 5*»3n0.* 
0'^“2)  f TC!H.(10,2)  = 3 

A (1  : , 3)  = ( . ItO  *.,+a  . '‘3  f35  .1  A + 3**2n*0 .2',»Sf3.  / (3.  »S-'*3»- 

ppino-’f^fjono)  s icoL(ia,3)=5 
AdO.'O  =-l.  3 ICOL  (Dt'tlrlO 


A(ll,l)  =23.37  7 ICCL ( 11, 1) =11 

4(11,-’)--  (1.  /3)  ' (P(l>  • ( .5*5*1 .)  *C?XP(  -.3*5)  ) 


3 100.(11,2) =20 


A(l?,l)  = .0'525  $ rCGL(l2,l)=4 

A(12,P)  =2.  5 ir0L(12,2)=12 

A (I  2 , ^)=  (!./'')*  P ( 1)  ♦ ( .5-S  H .)  *Ct;XP  ( - .3*3)  * ( ( .OOIO'S*  .0  174)  / ( .3* 
OSH.))  3i:OL  ( 12,3)=2C 

A (13  , 1)=  (1  .10  ..555) /r  ;*  IC  Cl.  ( 1 3 , 1 ) =2 

a(13,2)  = (.3'.»P)*3*.i.2*P()/(S*l.)  310  01  (13, 2)=3 

A(l3,3)  = ((.3'HA0-.O0  41)-S**2>(.42*AO-.10f.2)*S-.C0  21l/(3**2*3) 

•0  • ICO!  (13,  3)  =4 

A(13,O  = ((.143is*A:-3.*.17)*S**3*1.224  0*A0-l.riS)*S*’2*-(.U22*AQ- 
0.  .13'17)  *5  AO-.GIFf,)  / (4  .*3  ‘*3  )4  .♦?**  2*5)  $ 100'- ( 1 3 ,4)  =5 

4(13,*^)  = ( . 1 1 J7*5  * . C'5'53) /5  5 IC 01  ( 1 3, 5 ) = 1 2 

A(13,e,)  =(.33.j^‘3-.C3.,5)/S  $ ICOl  ( 1 3 , 6)  =13 

A(14,1)=-5TC  TTCOl (14,1)=1 

A(l'.  ,'>)=  VT1»0T:‘4-FC  • 1001  ( 14,2)  =2 

A(14,3)=  WO“OTa»OFO  7ICC1 ( 14, 3)  =3 


A(i; 


•»)=  -U. 


'.I  )*STC“  cro  'itcoi  (14,4)-:a 


A(i4,;)=  -JG‘C'o*SFj  ?rcti(i4,r.)=o 

A (14,!-,)=  (-U)  ‘S*;*vrc  *3  40*01:,  •CF.  *Gn*TFC)/TFC  310  0'.  (14  ,tS)  = 14 


Fig.  21.  (Continued)  Elements  of  the  A Matrix 


> , it 
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best  AVAllABU  COPY 


n( ' n , I ) =:caom  - rc  ‘OK'UMO’Tr 0**  j»r.c ) ( . c=’CTQ‘:?j »ni<  'K43» 

nr.-).*  ..-.C)  : M'.C'I.  (10, 1)  =1 

A (10  , > ) - '■RO'  ( • STO»s'‘G  -OS  j"cro)  •rTG»rc.o  ♦CTr»SF0‘^.T0*soa-.5* 

rrr- ^ OF fj'-*,-’)  i ictM.  ( lo  j.”) -2 

A ( ' ":  , Y)  = -F40  MO'::  *CTO  'FT  O'Or-J-CTO  •‘■'TO  ‘Sr-C*-.  O'FT  j'  0 ; « T FO  » * ? ) f . 

FTAT-M  crOO'S'C  ’ O'"!  SP  0 )•'•*(  J 'COO  fPTC'CFC  'PrC"SOC-.5'CT  0 ‘C^O  'OP* 

rr,'''T'^0"?)  f IFCL(!0,0)  = ? 

A(ir,',)=cAo*{-oRci-cTo»p'0'cr.j)*sA'j'((-sr>o'STa*cFfl»:>o*5Fij)'‘?ro* 
re  on  r tool  ( 1 0 ,4) 

A(io,::)=CA';»(-'',o'-sTo,*pio'or.j-cTO'RTo*-':r.op.o*rTC'n'>'oo»TF’.)*»?)f 
rOA  jMOSC'Cr  C'O'^u'PTG'CGC -ST  >J»CFO»XTC*OOC  *.0*510 'Crj'0“*3C-‘TKO-'»2) 
ir  'L  t rno)  =s 

A(1‘;  ,'■>)  :5A0M  (-r-onOTC*:  FCfSSnCPO)*PT)*cr,0-CTC*SriMT  J‘5Gf;  *.5* 
cr^:  * s-o*C)P*".r  * rr  3 $ icoi(15,f)=9 
Ad'.:,?')  ^(rAO'-src-SAO'CTO'CFj)  *(OP*r,(:'rF0-.5*o<'‘p.'n'03*vT0'r'0*‘2) 
TO  M.  (1  •},  n = l(. 

A(l';.,T)  =PPC»5LC  Z TOOL  (IF, 'Y)  = 15 

A { 15 , 1)  = ( -.p*  )'< ‘P'tO'TFC*  * T'OC)  *CA  0*CT  0 MFO  S ICO.  ( 15  , I ) = 1 
A (15,?)  --'^AC'  (?.Tr,  COO*  (('S'-dST  i'SP'o-rsc'TFn  MTO'SSO'oro'Sf'o-.o’Op* 
rr.r:T'^''o»?*rT'3<'SF0>  £ IC  OL  ( 1 5.  ? ) = 3 

A d ', , 3)  = - pAo  *■  ( 0 erd'  (5 sc'STO'S- O’-rsc-CFO) ) • i :?.  ( ir-,,  n = 7 

A do  ,.)  -0  a:  • (•J''.'*CG.d  (C5  0 'CTO' 5F0  ) -F  T . ‘ SGO  ‘ STC  *Sr  0 . 5 * 0 ” -GC  ^ T FJ  *»  2 
C*Pro*''n)  $ TPCL(16,5)  = d 

A ( l^^,'; ) = cc  0*  ( =>T  O'CGC*  ( Cl  C'ST  J'CFO  •■OrOMFO)  ♦■P.rO'SGO'  oro  *CF1-  .3*0P‘ 
roO*TFC",?'C’'0- CFO)  ■;  iC0Ld6,5)=9 

A d')  , j>  =rA:  * ( - OP  'GC'TFC*  P-TO'SFQ*  OK'PHO*  V rC*GC'r='C**2*CT0*5FQ) 

T''Tt.  (d  , 5)  =1 A 

AdO,?)  =-l.  I IC0Ld6,7)=16 


A d.:’ , 1 ) = -PPC*''  ? / ( ,?  . *7T  C**2)  -OP-*  *2 ‘T  FO  * » e/?  . * V T 0 * T r 0 ' * 2 * 0 P*0  < » -D'O* 
FTP)  5 1001.(17,1)  d 

Ad7,-n  =-T‘’»'2'7:o*-TFC*yT';'7FO**2*OP' OK*.PHO*VTQ  i I C CL  (1 7 , 2 ) "1 '• 
Ad',')  =-l.  £ ICOL  (17 ,3)  =17 


A d 3 , 1 ) = -rA  • ^ ?0  5.  * . 5095/7110 
A d ? , .’)  - C A ; ' 2 1 3 . ■ . 0 C 32  1 /7  TO  $ 

Ad  5 , ••)  - OAO'  ? i 3 .-'  ( . C 10  7 - . 021*10)  /'7TC 
A d 3 , .)  =CA';  * (-  ;c- FTC'S  F'.  ) /7T() 

Ad'),  -T)  =CAO»  ( TO'O'C'CFO)  /7T0  T 

A d 3 ,5)  =CAO  *2  iF  . • ] /'JTi 

Add,  /)  =rAi:  *73T.  ^ . 1312/7T0 
Ad3,n  =-i. 


$ ICOL  ( lf»,  1)  =2 
ICOL  (13,2)  ='» 

5 ICOL  (1.3, 3)  =& 
* ICOL (15, 5) =3 
ICOL(io,F) =9 
S rrCL ( 13,6)=12 
5 I'0L(18,7)=13 
■£  ICOL  (ie,a)=i3 


Ad9, 1)  =PPO*CLO 
Ad9,2)  =-l. 


« ICOL (10, 1) =15 

J IOOL( 19,2) =19 


Af'»3,l)  ='<0/»3'’.1  5 ICOL(20,1)  = 16 

A(2J  , 2)  =-'<!; ''■"T'TP>AO/PD?U  * ICOL  ( 20 , 2)  =17 

A('0,n  =-<3  • )r!l)/PP'2C  ■£■  ICOl.(  20,3)  =ia 

A(?0,'O  r-i':  • (O'Md-nrnC' Hn'tFGAJ) /R(120"2  $ rOJL<2o  ,'0  =19 

A(20,'i)  =-l.  3 irOL(20,';)=20 


Fig.  21.  (Continued)  Elements  of  the  A Matrix 
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h-  5C  cNi 

JCCQQ:'*-  C1.^-3Q<J 
3aaacLcrJ--9-a3-o->o'0‘0‘o+j^uxJ  jcxcl  , 


* -X  * -K 


« « X « « « 


XXX-X*-tt*-»i 


X * X ■X 


X X X X 


X X X X 
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BEST-  AVAIUBLE  COPY 


n(t)  (1.«0*Sft0-.8O8»A0*CA0+.0O0  5*Cft0-.0O05M0»SA0f2.  325*A0 

r»»'»»SAO)  ) /VTO)  • 9(8) 

n( ’)---((  /VTO-(3.9’-6-8.204»AO)  • ((-S/VT0)/(l.*(3.*f:»‘5)/(?I* 

rv*l>  ))>  •“('») 

Tf  •')  =-  ( * (- J.  80»CA0-.  36  3'A8*SAOf  . 00  05  • AO 'CA  0 «■ . 89  05  *S  A 0-2 . 325* 

rA0**2»CA0) ) /VTO) 'PCS) 

9(..)  = ( ( »1  .?9*3.,TT5*A0)  /VTO  - ( . 020<.*2.033»AC)*  ( (-5/vrO)  / (1  .♦<3.*E 
r*<^)/<°I»VT0))))*P('‘>*-.l02''P(6) 

'>(5)  = ((»i.(,7)/VTO-(-.2?l*.0/ir)*((S/VTO)/(l.*(4.»''*S)/<PI*VrO 
f”)  ) ) ) »“  <5) 

P('))  = -{(2.<.3-r,  .l'TAO)/VTO  -(-.  1063)*(  (-S/VT0)/(1.*(3.*E*S)/{PI* 
r.V’O  )))) 'P  ('*)♦.  0 12/8»PL6) 

! 1 > = f.  0C?''»5»  2*.  0 526*5*^  .0100)  / (':•*  3V5.*5**2)  *P(  2) 

r*'((-3.'>85*(-3.10'’‘CA0-.?f.".»A0*3A0».090  5‘A}*CA0f.0O0  5»SA0-2.32  5*AO 
C*»•’»C^O)*(.G/8fi’5♦*3*■.3"^T^*S**^t.830*S^•.2■'?8)/(S♦'3^30.*S♦♦2)) 
r;/VT0)*n(5) 

P ( 1 1 ) = t . • P ( 3 ) 

p(12)  = (.001'3*5v.0  075),/(.3*-SH.)*P(3) 

n{13)  = (((-(l.ll’St^.553)/S)/7TO)  -((.O13-.03/*AO)*5«-(.nO7-.019*AO)) 

r,/5)*((-S/VT0)/(l.K3.»E‘5)/(PI*V10))))*P('*)*(  ( . 0 0'*»  3 ♦ . 0 0 2 ) /S)  » P (6) 


Fig,  23.  Elements  of  the  B Matrix  (Basic  F-15) 


n(t ) =-  ((286  .*  (3 .30*SA0-.  863*A0*CA0*  . 0905*CA0-.0O05»  AO'SA  04-2. 325*A0 
C**2*SA0) ) /VTO) *P(5) 

P(2)i)=  ( ( ♦245.8)  / VTO- (2 .94  6 -3. 29'**  AO)  M(-S/VT0)/(l.f(3.*E*S)/(PI» 
CVTO) ) ) )*P(4)  -.913*0(6) 

P(3) =- ( (386.* (-3.  30*CA0- .b68*A0*SA0* .0905*A0*CAn* . 0905  »SAO-2,325* 
rA0**2*CA0)  ) /VTO  ) ’‘P(5) 

9(4) =( ( H .29  *3. 3 35*A0) /VTO  - (.0204*2.833*A0)* ( (-S/VTC) /(I . ♦(3.*e 
C*S) /(PI*VTO) )) ) *P(4) ♦.1C2*P(6) 

9(5)  -(  (♦1.47)  / VTO  - (-.221  ♦ . C717)*{  (S/'/TO)  / (I.  *(4  . *P’ 3 ) / ( PI  * VTO 
C);.> ) ) * P ( 5 ) 

9(6)= -((2.43 -6. lO'AO) /VTO -(-.1063 )♦( (-S/VrO) /(!.♦( 5.*E»S)/(Pt* 
CV^O) ) ) ) *P(4) ♦.01276*P(6) 

'9(lO)=((-8.335*(-3.802*CAO-.368*AO*3AO^. CP05*A0*nA6 ♦ .0905*SAO 
C-r’.  325*A0•*^*CA0)*(.0166*S•*2♦.C4  9*S♦.^3  55)/(30.‘S*♦^^S**3))/VT0 
0*0(5) 

^ 9(ll)=l.*P(3> 

0(12)= (.0019*S*.0074)/(.3*S*1.) *P(3) 

9(l3)=(((-(1.03*?^.52)/9)/VTO*((.012-.C35*AQ)*3*(.OC6-.Ol7*AO))/S 

O*((-S/VT0)/(l.*(!.*E»S)/(PT*VTC))))*=>(4)*{(.0030*Sv.0019)/S)*P(r>) 


Fig.  24.  Elements  of  the  B Matrix  (Automatic  Fire  Control  System) 
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